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ABSTRACT 

Various  theories  of  physical  adsorption  are 
reviewed.  Although  some  theories  which  treat  mobile  adsorp¬ 
tion.  are  discussed,  the  emphasis  is  on  localized  adsorption 
theories.  The  work  of  Ballentine  and  Betts  in  which  the 
Ising  model  of  ferromagnetism  is  applied  to  the  problem  of 
localized  monolayer  adsorption  with  attractive  lateral 
interactions  is  examined  in  some  detail.  An  extension  of 
this  in  which  Ballentine  applies  the  three  dimensional  Ising 
model  to  multilayer  adsorption  is  outlined,  and  results 
which  show  the  variation  in  the  location  of  the  critical 
point  with  the  number  of  layers  are  given. 

A  model  for  two  layer  adsorption  with  attractive 
lateral  interactions  is  proposed.  This  model  is  suggested 
by  the  Ising  model,  but  it  is  more  realistic  concerning  the 
nature  of  the  second  layer  sites.  It  is  essentially  a  two 
dimensional  model  in  which  each  adsorption  site  is  considered 
to  be  intone  of  three  possible  "occupation  states.”  The 
model  is  specialized  in  order  to  simplify  the  mathematical 
treatment.  The  partition  function  for  the  system  is 
expanded  as  an  exact  series  in  both  high  and  low  temperature 
variables.  The  calculation  of  the  series  coefficients 


- 
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involves  the  technique  of  configuration  counting.  An  estimate 
of  the  location  of  the  critical  point  for  this  system  is  made. 

The  work  of  the  thesis  is  summarized  and  the 
possibilities  of  generalizing  the  treatment  and  extending 
the  model  to  include  more  layers  are  discussed. 
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Chapter  I.  Introduction 


1 . 1  Adsorption  -  Definition  and  History 

When  a  gas  or  vapour  comes  into  contact  with  a 
solid,  the  gas  may  be  adsorbed  onto  the  surface  of  the  solid 
or  absorbed  into  the  interior  of  the  solid.  Regardless  of 
the  type  of  force  or  combination  of  forces,  a  molecule 
of  gas  absorbed  into  a  solid  is  subject  to  forces  equally 
in  all  directions.  On  the  other  hand  molecules  adsorbed 
onto  the  surface  are  subject  to  unbalanced  forces  which  hold 
it  onto  the  surface.  In  effect  the  solid  has  a  "surface 
tension"  and  the  presence  of  the  adsorbed  molecule  reduces 
the  surface  area  and  lowers  the  surface  free  energy.  Thus 
adsorption  is  an  exothermic  reaction  and  the  so-called  heat 
of  adsorption  is  a  physically  measurable  quantity. 

Both  absorption  and  adsorption  can  take  place 
through  either  a  chemical  or  physical  process.  If  there  is 
a  transfer  of  electrons  between  the  atoms  of  the  gas  and 
the  solid  then  a  new  compound  is  formed  the  process  is  said 
to  be  chemical  ab-  or  adsorption.  If  however  the  forces 
are  of  the  Van  der  Waals  type,  i.e.,  strictly  molecular 
interaction,  then  the  process  is  physical.  In  physical 
absorption  a  solid  solution  is  formed.  In  the  process  of 
physical  adsorption  the  adsorbed  film  behaves  in  many 
respects  like  a  two-dimensional  liquid.  Physical  adsorption 
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may  further  be  subdivided  into  localized  and  mobile  adsorp¬ 
tion.  In  the  former  the  gas  molecules  which  adhere  to  the 
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solid  surface  are  essentially  bound  to  fixed  adsorption 
centers ,  whereas  in  the  latter  a  thin  liquid  film  is  formed 
and  the  molecules  are  free  to  move  within  the  film.  The 
model  which  we  will  propose  later  in  this  thesis  will  be 
for  localized  physical  adsorption  with  lateral  interactions. 

Adsorption  was  first  discovered  in  1773  by  Scheele 
who  found  that  charcoal  in  contact  with  a  gas  in  some  way 
reduced  the  amount  of  gas  in  the  free  state.  The  first 
person  to  perform  a  series  of  systematic  experiments  on 
adsorption  was  de  Saussure  (l)  who,  in  l8l4,  measured  the 
amount  of  adsorption  of  various  gases  on  many  different 
adsorbents.  Since  then  large  numbers  of  experiments  have 
measured  the  adsorption  of  nearly  every  conceivable  combina¬ 
tion  of  adsorbate  (gas)  and  adsorbent  (solid).  A  great 
difficulty  in  the  study  of  adsorption  has  been  that  many 
conflicting  results  have  been  reported.  The  reason  for  this 
is  that  the  physical  properties  of  the  adsorbent  depend  so 
much  on  the  way  in  which  the  sample  is  prepared.  Thus  samples 
of  the  same  adsorbent  prepared  in  different  laboratories, 
and  sometimes  even  two  samples  from  the  same  laboratory,  may 
give  far  from  similar  results.  It  was  pointed  out  later  that 
many  of  the  earlier  experiments  were  performed  using  samples 
of  adsorbent  which  already  contained  a  composite  layer  of 
adsorbed  gas  molecules  and  had  not  been  properly  "cleaned".  This 
effectively  suppressed  the  properties  of  the  intended  adsorbent. 
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These  factors  made  it  difficult  to  propose  and 
check  a  theory  of  adsorption  since  there  were  few  consistent 
experimental  results.  However  in  1911  and  1915  Knudsen  (2) 
showed  that  if  one  assumed  the  dimensions  of  the  apparatus 
were  sufficiently  small  with  respect  to  the  mean  free  path 
of  the  gas  molecules,  one  could  do  calculations  using  the 
kinetic  theory  of  gases.  Knudsen Ts  experiments  tended  to 
support  his  theoretical  calculations.  Since  that  time 
various  theories  have  been  proposed  and  we  shall  discuss 
them  later  in  this  thesis. 

1 . 2  Cooperative  Phenomena 

We  consider  the  class  of  statistical  mechanical 
systems  whose  macroscopic  properties  are  determined  primarily 
by  the  interactions  between  the  microscopic  subsystems. 

Systems  in  this  class  are  said  to  exhibit  cooperative  phenomena. 
Other  than  the  Bose  gas,  all  systems  which  show  phase  transi¬ 
tions  fall  into  this  class.  The  physical  adsorption  of 
molecules  which  interact  laterally  is  a  cooperative  phenomenon. 
This  exhibits  a  phase  transition  known  as  critical  adsorption, 
where  below  a  given  temperature  the  number  of  adsorbed 
molecules  increases  sharply  at  a  particular  pressure. 

The  simplest  model  of  cooperative  behaviour  was 
proposed  in  1925  by  Ising  (3).  In  this  model  Ising  assumed 
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that  all  subsystems  of  a  cooperative  assembly  were  identical 
and  could  be  in  one  of  two  possible  states.  He  also  assumed 
that  the  interaction  between  subsystems  was  sufficiently 
short  range  that  all  but  nearest  neighbour  interactions 
could  be  neglected.  Ising  applied  this  model  to  the  problem 
of  ferromagnetism.  Since  then  other  researchers  have  extended 
the  theory  and  applied  it  to  other  cooperative  systems. 

Fowler  (4)  and  Peierls  (5)  were  the  first  (1936)  to  apply 
the  two-dimensional  Ising  model  to  monolayer  adsorption  with 
lateral  interactions.  However,  their  solutions  made  use 
of  two  different  approximations  and,  although  they  obtained 
a  critical  point,  their  results  in  this  region  were  very 
distorted. 

As  recently  as  i960  Bumble  and  Honig  (6)  used 
another  approximate  method  (due  to  Hijmans  and  de  Boer  (7)  ) 
to  solve  this  same  problem.  In  1962  Ballentine  and  Betts 
(8,9)  attacked  the  problem  using  exact  series  expansions 
and  showed  that  complete  adsorption  isotherms  could  be 
calculated  in  this  way.  Ballentine  extended  the  approach  to 
include  two  or  more  layers  of  an  Ising  model  and  found 
estimates  of  the  critical  adsorption  temperature. 

1 . 3  Reasons  for  Studying  the  Adsorption  Problem 

Ballentine  in  the  introduction  to  his  thesis  has 
outlined  the  reasons  for  which  this  research  was  started. 
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Various  experimentalists  have  reported,  adsorption  isotherms 
which  exhibit  a  critical  point.  Jura  et  al  (10),  Fisher  and 
McMillan  (ll),  and  Machin  and  Ross  (12)  have  all  found  iso¬ 
therms  with  single  di scontinuities ,  whereas  Singleton  and 
Halsey  (13)  have  found  a  stepwise  isotherm  which  has  two 
distinct  discontinuities.  This  latter  result  is  for  the 
system  krypton  on  graphitized  carbon  black. 

The  work  done  by  Ballentine  was  very  successful 
in  providing  isotherms  for  critical  monolayer  adsorption. 
Although  his  study  of  the  n-layer  adsorption  problem  by 
direct  application  of  the  Ising  model  is  certainly  justified 
for  the  reasons  he  states,  the  model  does  not  seem  physically 
realistic  in  regard  to  the  nature  of  the  sites  for  adsorption 
in  the  second  and  higher  layers.  The  model  which  we  will 
propose  here,  although  exploited  hardly  at  all,  will  be 
an  improvement  in  this  regard.  The  Ising  model,  by  its  very 
nature,  will  show  but  one  critical  point  regardless  of  the 
number  of  adsorbed  layers.  Physically,  however,  for  an 
ideally  homogeneous  adsorbent  one  would  expect  critical 
adsorption  as  each  layer  is  filled.  This  is  the  way  in 
which  Halsey  and  Singleton  interpret  their  result,  and  this 
is  also  supported  by  the  work  of  Polley,  Schaeffer,  and 
Smith  (l4)  who  showed  that  the  larger  the  crystals  in  the 
adsorbent  the  more  step-like  the  isotherm  becomes.  A  further 
reason  for  proposing  this  model  is  that  if  it  can  be  solved 
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without  the  restrictions  which  have  been  placed  on  the 
interaction  energies.,  then  it  is  believed  that  stepwise 
isotherms  would  result  --  each  layer  filling  discontinuously 
as  the  pressure  is  increased  or  the  temperature  lowered. 

1.4  The  Contents  of  the  Thesis 

In  chapter  II  of  this  thesis  we  will  survey  many 
of  the  early  and  current  theories  of  adsorption,,  starting 
with  the  rather  simple  approach  of  Henry Ts  law  and  Langmuir Ts 
isotherm,  and  continuing  on  to  the  popular  theory  due  to 
Brunauer,  Emmett,  and  Teller  and  some  of  the  modifications 
that  have  been  suggested  for  it.  Chapter  III  will  contain 
an  outline  of  the  work  done  by  Ballentine  and  Betts  for  both 
monolayer  and  n-layer  adsorption.  A  reasonable  understanding 
of  their  work  is  recommended,  for  it  was  this  which  suggested 
the  model  which  is  proposed  in  Chapter  IV.  We  will  of 
necessity  specialize  the  model  and  from  the  series  expansions 
for  the  partition  function  estimate  the  critical  temperature 
for  the  adsorbed  system.  In  the  concluding  chapter  we  will 
compare  this  result  with  those  obtained  from  other  theories, 
and  also  discuss  the  possibilities  which  remain  for  the  more 
general  model. 
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Chapter  2.  Other  Theories  of  Adsorption 


2 . 1  Monolayer  Adsorption  Theories 

Among  the  many  theories  that  have  been  proposed 
to  explain  physical  adsorption,  the  simplest  is  HenryTs 
Law.  It  assumes  that  the  amount  of  gas  adsorbed  is  directly 
proportional  to  the  equilibrium  pressure  of  the  gas.  To 
derive  this  law  theoretically  one  must  assume  that  the 
behaviour  of  the  gaseous  and  adsorbed  phases  is  perfect. 

This  is  possible  only  if  both  phases  are  extremely  dilute. 
The  linear  isotherm  which  describes  the  adsorption  can  be 
written  as 


v 


ad  2  e  /  RT 

p  rt  sq  e 


(2.1) 


where  e  is  the  adsorption  potential,  a  is  the  fraction 
of  the  adsorbent  available,  d  is  the  thickness  of  the  film, 
o>  and  are  in  the  same  ratio  as  the  internal  energies  of 

molecules  in  the  gas  and  adsorbed  phases  respectively,  and 
p  and  v  are  the  equilibrium  pressure  and  the  volume  adsorbed 
respectively.  Isotherms  which  appear  to  be  linear  have  been 
found  experimentally,  but  they  cover  only  a  very  small  part 
of  the  possible  pressure  range.  Since  the  pressure  range 
is  so  small,  and  since  more  accurate  measurements  on  these 
same  systems  now  show  slight  deviations  from  linearity,  it 
is  widely  held  that  this  picture  is  far  too  simple  to  be  a 
valid  theory  for  adsorption. 
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The  next  theory.,  in  order  of  ascending  sophisti¬ 
cation,,  is  that  due  to  Langmuir  (15) •  Although  nothing  is 
said  about  the  mechanism  of  adsorption  in  Henry’s  Law,  we 
are  able  to  see  that  the  adsorbed  layer  is  mobile,  there 
being  nothing  to  hold  a  particular  molecule  to  a  particular 
spot  on  the  surface.  On  the  other  hand  in  Langmuir’s  theory 
we  assume  that  the  gas  molecules  are  adsorbed  onto  the 
solid  at  definite  sites,  and  that  each  site  is  capable  of 
holding  only  one  gas  molecule.  We  must  also  assume  that 
there  is  no  interaction  between  the  adsorbed  molecules. 
Originally  the  isotherm  equation  was  derived  kinetically  by 
assuming  that  the  gas  was  in  equilibrium  with  the  adsorbed 
layer.  By  equating  the  rate  at  which  the  molecules  evaporated 
from  the  surface  to  the  rate  at  which  they  "condensed" 
onto  the  surface,  Langmuir  in  1918  obtained  an  isotherm  of 
the  form 


e 


bp 

1  +  bp 


(2.2) 


where  6  is  the  fractional  coverage,  p  the  equilibrium 
pressure,  and  b  a  function  only  of  temperature.  In  1935 
Fowler  (16)  gave  a  statistical  derivation  of  equation  (2.2). 
It  is  interesting  to  note  that  for  very  low  pressures  (i.e., 
when  both  the  gaseous  and  adsorbed  phases  are  very  dilute) 
Langmuir’s  isotherm  reduces  to  Henry’s  Law.  For  the  other 
limiting  case  when  p  00 ,  we  find  Q  -»  1  or  v  -*•  vmax^  where 
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v  equals  the  volume  of  a  completely  filled  monolayer. 

maA 

This  model  does  not  permit  any  adsorption  beyond  a  monolayer. 

Williams  ( 1 7 )  and  Henry  (l8)  have  modified  the 
Langmuir  model  slightly  by  considering  the  case  where  a 
gas  molecule  when  adsorbed  occupies  n  adjacent  adsorption 
sites  on  the  lattice.  Several  others  have  derived  an  isotherm 
for  the  case  where  the  adsorbed  molecule  dissociates  into 
n  atoms  and  each  atom  occupies  a  site.  However  these 
isotherms  are  not  very  important  in  physical  adsorption 
and  won’t  be  considered  further. 

Volmer  (19)  has  treated  the  case  for  mobile  ad¬ 
sorption,  i . e . ,  where  the  adsorbed  molecules  are  free  to 
move  in  two  dimensions  on  the  surface  of  the  adsorbent. 

By  assuming  that  the  gas  was  perfect  and  that  the  adsorbed 
phase  was  a  two-dimensional  gas  obeying  the  (non-perfect) 
equation  of  state 

b(A  -  N  £)  =  NkTj 

where  b  is  the  surface  pressure,  A  the  area,  and  £ 
a  constant,  Volmer  obtained  an  isotherm  of  the  form 


N 

A 


_ 

1  +  b!p 


(2.3) 
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which  is  similar  in  form,  but  not  identical,  to  the  Langmuir 
isotherm  (2.2)  for  localized  adsorption. 

Fowler  (4)  in  1936  made  a  very  important  contribu¬ 
tion  to  the  theory  of  adsorption  when  he  proposed  a  model 
which  took  into  account  the  attractive  lateral  interactions 
between  molecules  in  the  adsorbed  phase.  His  model  was 
simply  the  Ising  model  of  ferromagnetism  applied  to  adsorption. 
He  said  that  the  adsorption  was  localized  and  the  attraction 
was  constant  and  sufficiently  short  range  in  nature  that 
all  but  nearest  neighbour  interactions  might  be  neglected. 

This  model  will  be  dealt  with  in  much  greater  detail  in 
chapter  III. 


Since  at  that  time  no  exact  methods  of  solution 
were  known  for  the  Ising  model,  Fowler  made  use  of  the 
Bragg -Williams  (20)  approximation  in  his  solution.  The 
number  of  nearest  neighbour  pairs  of  adsorbed  molecules, 
which  is  required  in  the  calculation,  is  in  this  approxima¬ 
tion  taken  as  the  average  number  that  would  occur  for  a 
given  coverage  if  the  distribution  were  random.  On  this 
basis  Fowler  derived  the  following  isotherm  equation 


bp 


e 


-e  V/  KT 

•> 


(2.4) 


where  0,  p,  and  b  have  their  usual  meaning,  -V  is  the 
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constant  interaction  energy  between  nearest  neighbours,,  and 
q  is  the  coordination  number  (or  number  of  n-n  pairs  per  site) 
for  the  lattice.  The  most  noteworthy  property  of  this 
solution  is  that  it  exhibits  critical  adsorption.  The  speci¬ 
fic  heat  is  discontinuous  as  a  function  of  temperature  and 
the  isotherms  for  temperatures  below  a  critical  temperature 

T  have  discontinuities.  For  the  Bragg -Williams  approx! - 
c 

mation  the  critical  temperature  is  given  by 


V 

KTC 


1  , 


(2.5) 


for  a  simple  quadratic  lattice.  If  one  applies  the  quasi¬ 
chemical  or  Bethe  approximation  to  the  problem  with  lateral 
interactions  (5)  one  gets  a  more  accurate  value  of  the 
critical  temperature 

U-  =  2  In  2  =  1.3863  •  (2.6) 

1 C 

We  shall  see  in  the  next  chapter  how  the  critical  temperature 
has  been  found  exactly. 


2.2  Two  Layer  Adsorption  with  Critical  Conditions 

The  work  of  Cernuschi  (21 )  and  Dube  (22)  is  cer¬ 
tainly  worthy  of  mention  in  this  review  of  earlier  adsorption 
theories.  Each  of  them  used  a  two-layer  model  in  which 
four  adsorbed  molecules  on  adjacent  sites  (forming  a  square) 
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in  the  first  layer  constitute  a  site  for  the  possible 
adsorption  of  a  molecule  onto  the  second  layer.  Cernuschi 
neglected  lateral  interactions  while  Dube  included  them. 

Both  of  them  found  results  which  indicated  that  critical 
adsorption  could  take  place  for  their  models,  providing 
the  binding  energies  satisfied  a  certain  condition. 

Cernuschi  considered  the  molecules  in  the  first 
layer  to  be  bound  to  the  surface  by  an  energy  -  xo  *  and 
since  there  was  no  interaction  all  configurations  were 
equally  probable.  On  this  basis  the  average  number  of 
available  sites  in  the  second  layer  is  (N./M^)  ,  where 

=  the  number  of  molecules  in  the  first  layer,  and 
=  the  number  of  sites  in  the  first  layer.  The  molecules 
in  the  second  layer  are  bound  to  the  surface  with  energy 
-  .  After  finding  the  partition  function  using  these 

assumptions,  Cernuschi  showed  that  critical  adsorption  takes 
place  subject  to  the  conditions 

Sn  =  i  ,  and  KT  =  — - -  .  (2.7) 

lc  4  c  0#42 

Thus  if  y  >  6  no  critical  point  exists.  is  the 

^o  o  lc 

critical  coverage  of  the  first  layer. 

Cernuschi  gave  a  brief  discussion  of  part  of  the 
effect  of  an  attractive  interaction  x^  on  this  model. 

He  neglected  the  effect  of  the  interaction  on  the  existence 


♦ 
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of  second  layer  sites,  and  used  the  Bragg -Williams  approxi¬ 
mation  to  calculate  the  increased  binding  energy  in  the 
first  layer  due  to  the  interaction.  Cernuschi^  results 
showed  that  the  first  layer  critical  coverage  decreased 
with  increasing  attractive  interaction  between  molecules, 
while  the  critical  temperature  (for  fixed  xo  and  6  ) 
decreased  to  a  minimum  for  some  value  of  the  interaction 
'(Xq  ~  3(§0  -  XQ )  )  and  then  increased  with  increasing  Xq • 

Dube  used  exactly  the  same  picture  of  adsorption 

but  did  better  calculations  using  the  attractive  interaction 

Xq •  He  follows  Peierls  (5)  treatment  using  the  Bethe 

approximation,  rather  than  the  Bragg -Williams  method  as 

Fowler  did,  to  calculate  the  effect  of  the  interaction  on 

increasing  the  probability  of  existence  of  second  layer  sites 

The  equations  he  derived  had  to  be  solved  numerically  for 

different  values  of  the  interaction.  He  defined  the 

probability  of  existence  of  a  second  layer  site  as  0^Mq  . 

Cernuschi  used  0^M1  =  (N1/M1)^M1  for  this  quantity.  Dube 

showed  that  for  all  values  of  Xq  corresponding  to  a  non- 

4 

zero  attractive  interaction  0^  >  0-^  ,  except  at  the  end 

4 

points  @1  =  0  and  0^-1  where  0^  =  @1  .  Dube  then 
used  the  Bragg -Williams  approximation  to  set  up  the  partition 
function  for  the  two-layer  adsorbed  system.  He  was  then 
able  to  find  all  the  critical  functions  in  terms  of  the 


14. 


energies  6q,  xqj  X j_’  an<^  a  parameter  p  which  enters  in 

the  functional  fit  of  0^(0^,  x^)  ^o  the  numerical  results. 

DubeTs  results  again  showed  that  the  necessary  condition 

for  critical  behaviour  was  that  5  >  y  . 

o  ^o 

PowlerTs  results,  using  the  Ising  model,  for 
monolayer  adsorption  show  that  the  presence  of  attractive 
interactions  between  adsorbed  molecules  is  sufficient  to 
cause  critical  adsorption  to  occur;  while  according  to  the 
model  of  Cernuschi  and  Dube  critical  adsorption  will  occur 
in  two  layer  adsorption,  whether  there  is  lateral  interaction 
between  molecules  or  not,  so  long  as  the  second  layer  molecules 
are  bound  more  tightly  to  the  surface  than  those  in  the 
first  layer. 

2 . 3  Multilayer  Adsorption  Theories 

Much  of  the  early  theoretical  work  in  adsorption 
treated  one  or  possibly  two  layers  of  adsorbed  molecules, 
rather  than  many  layers.  The  reason  for  this  was  that  most 
early  experimental  work  was  done  on  charcoal,  which  is  a 
unique  adsorbent  in  that  adsorption  onto  it  rarely  occurs 
beyond  a  monolayer. 

The  first  theory  of  multilayer  adsorption,  known 
as  the  potential  theory,  was  introduced  in  1914  by  Polanyi 
(23).  This  theory  is  still  accepted  today,  largely  because 
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it  does  not  insist  on  a  detailed  picture  of  the  physical 
process  involved.  It  treats  a  mobile  film  by  considering 
it  to  be  trapped  in  the  potential  gradient  of  the  surface 
force  which  decays  with  the  distance  away  from  the  surface. 

The  adsorption  potential  e  is  the  measure  of  the  force 
on  an  adsorbed  molecule,,  the  adsorption  potential  being  the 
work  required  to  bring  a  molecule  from  the  ground  state  of 
the  gas  phase  to  the  point  in  the  film  at  which  it  is  measured. 
We  denote  the  volume  of  the  adsorbed  film  as  d.  We  assume 
that  e  is  independent  of  T  and  thus  that  the  curve 
e  =  f(d)  is  sufficient  to  describe  the  buildup  of  the  film. 

€  =  f(d)  is  in  essence  a  potential  distribution  function 
and  it  is  known  as  the  characteristic  curve.  This  curve 
cannot  be  measured  directly  from  any  experimental  system. 
Rather,  one  must  take  a  "master  isotherm",  one  which  is 
known  over  a  large  range  of  adsorbed  volumes,  and  calculate 
f(d),  and  then  from  this  curve  calculate  any  other  isotherm 
for  the  system. 

An  interesting  point  about  this  theory  is  that  for 
different  gas  systems  being  adsorbed  onto  the  same  solid 
one  can  make  all  the  characteristic  curves  coincide  with  one 
another  by  writing 


£  =  P  f(<t>)  • 


(2.8) 
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The  factor  £  is  known  as  the  "coefficient  of  affinity" 
and  is  different  for  each  adsorbate.  It  has  been  found  that 
the  physical  property  of  all  the  adsorbates  which  is  most 
closely  porportional  to  (3  is  the  molar  volume  of  the 
adsorbate  in  the  liquid  phase.  With  this  knowledge  one  is 
able  to  predict  something  of  the  adsorbability  of  a  substance 
without  actually  performing  any  experiments. 

In  1938  Brunauer,  Emmett,  and  Teller  (24)  success¬ 
fully  extended  the  Langmuir  model  of  localized  monolayer 
adsorption  to  solve  the  problem  of  multilayer  adsorption. 

The  BET  theory,  as  it  is  known,  made  certain  simplifying 
assumptions  that  made  it  possible  to  treat  the  problem 
mathematically.  First  it  is  assumed  that  molecules  are 
adsorbed  vertically,  a  molecule  in  a  given  layer  serving  as 
a  site  for  a  molecule  to  be  adsorbed  into  the  next  layer. 

All  interactions  between  molecules  in  the  film  are  neglected. 
One  denotes  the  areas  of  the  surface  covered  by  exactly 
0,  1,  2,  3*  •  •  •  i;  layers  of  molecules  by  sQ,  s^,  s^,  s^ 

.  .  .  si,  and  assumes  that  at  equilibrium  the  rate  of 
condensation  on  top  of  the  (i  -  l)—  layer  to  form  i  layers 
is  equal  to  the  rate  at  which  molecules  evaporate  from  the  i— 
layer.  The  equations  which  represent  this  equilibrium  are: 
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-E  /RT 

al  P  s0  =  bl  S1  e 

-Eg/RT 

a2  P  S1  =  b2  s2  e 


a. 

1 


P  s 


i-1 


(2.9) 


where  p  is  the  pressure  at  equilibrium,  E^ 

t  h 

adsorption  in  the  i —  layer,  and  a^  and  b_^ 


is  the  heat  of 
are  constants. 


The  total  surface  area  of  the  adsorbed  system  is 


given  by 


00 


A  =  2 

i=0 


(2.10) 


and  the  total  volume  adsorbed  is  given  by 


v 


v 


o 


00 

2  i 
1=0 


(2.11) 


where  vQ  =  the  volume/unit  area/monolayer. 

One  can  then  write  the  coverage  in  terms  of  the  number  of 
layers  by  00 


v 


v 


v 


m 


A  v 


o 


2  is. 

i=0  1 

00 

2  S. 
i=0  1 


(2.12) 


where  v  is  the  volume  of  one  complete  monolayer, 
m 
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The  assumptions  made  by  Brunauer,  Emmett  and 
Teller  which  effected  the  solution  were  that 


and 


=  E^  =  E^  =  . . 


b2  _  b3 

a^  a0  a 


3 


4 


=  E.  =  Et 
l  L 


b 


(2.13) 


This  is  effectively  saying  that  the  properties  of  the  second 
and  higher  layers  with  respect  to  evaporation  and  condensa¬ 
tion  are  identical  to  those  properties  of  the  bulk  liquid, 
for  by  E^  we  refer  to  the  heat  of  liquefaction  of  the 
adsorbate.  If  we  use  (2.13)  and  the  definitions 


x  = 


a. 

i 


b 


E/RT 

p  e  ,  i  /  1 


i 


y  = 


a  E/RT 

—  P  e 
1 


(2.14) 


we 


can  write  the  equations  (2.9)  as 


sr  y  so 


s2  =  X  Si  =  c  X  s0 


si  =  x  si-l  =  C  X  S0 


(2.91) 


where  c  =  y/x  .  When  these  are  substituted  into  the  series 
in  (2.12)  one  is  able  to  simplify  the  expression,  the 
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denominator  containing  the  sum  of  the  infinite  geometric 

series  §  x  ,  and  the  numerator  containing  the  product  of 

i=l 

x  and  the  derivative  of  that  series.  Equation  (2.12) 
becomes 

P  v  _  cx  , 

6  -  vm  -  Ti-x) (1  -  x  +  cx)  (2-] 

One  can  see  from  the  denominator  of  (2.15)  that  to  satisfy 
the  condition  that  the  volume  adsorbed  at  the  saturation 
pressure  be  infinite,  x  must  equal  one.  One  has  then 


(2.1 6) 


in  general,  and  the  substitution  of  (2.l6)  into  (2.15)  gives 


0  ”  v  “  (P  -  P)  il  +  ( c -1 )  P/P  ) 
m  v  o  '  o 

which  is  the  final  form  of  the  BET  "  00  -  form"  isotherm. 

Brunauer,  Emmett,  and  Teller  have  also  considered 
the  case  where  adsorption,  for  some  reason,  is  not  possible 
beyond  the  n—  layer,  rather  than  the  possibility  of  an 
infinite  number  of  layers  as  before.  For  this  case  the 
series  in  equation  (2.10),  (2.1l)  and  thus  (2.12)  must  be 

cut  off  at  n.  When  this  is  done,  rather  than  (2.15), 
one  gets  the  following  isotherm  equation 
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(9  _  v _ (  cx  x  /I  -  (n-l)  xn  +  nxn+^ 

Vm  1  '  X  1  +  (c-1)  x  -  cxn+l)*  f2-1?) 

which  is  known  as  the  n-layer  BET  isotherm.  It  is  interesting 
to  note  that  if  n  =  1,  (2.17)  reduces  to  the  Langmuir 
monolayer  isotherm.,  and  if  n  ->  (2.17)  becomes  the  same 

as  (2.15)*  both  of  which  results  one  would  expect  for  this 
theory. 

The  BET  theory  has  had  very  good  success  both  in 
fitting  experimental  isotherms,,  and  in  calculating  the  heat 
of  adsorption  and  the  total  surface  area  of  an  adsorbent 
from  experimental  isotherms.  However  the  theory  does  have 
its  shortcomings  -  mainly  its  neglect  of  lateral  attraction 
between  adsorbed  molecules,  and  its  assumption  that  molecules 
are  adsorbed  vertically  on  top  of  one  another.  We  shall  see 
that  various  workers  have  modified  the  theory  in  various 
ways  to  at  least  partially  overcome  these  objections. 

Halsey  (25)  has  considered  a  modification  in  which 
second  layer  molecules  are  adsorbed  onto  the  centre  of  a 
triangle  of  molecules  in  the  first  layer.  His  calculations 
show  that  unless  the  first  layer  is  over  l/3  filled,  there 
will  be  very  little  filling  of  the  second  layer;  and  in  any 
case  there  will  be  very  little  filling  above  the  second  layer 
unless  the  pressure  is  very  close  to  saturation.  We  will 
consider  HillTs  modification,  in  which  he  includes  lateral 
interactions,  in  the  next  section. 


2.4  Multilayer  Adsorption  with  Stepwise  Isotherms 

In  1947.,  Hill  (2 6)  gave  a  modified  treatment  of 
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the  BET  theory  in  which  he  assumed  that  lateral  interactions 

were  present  in  all  layers  but  were  of  sufficiently  short 

range  that  those  between  nearest  neighbours  only  need  be 

considered.  His  work  was  actually  to  the  BET  theory,  what 

Fowler Ts  was  to  the  Langmuir  theory,  for  Hill  used  the 

Bragg -Williams  approximation  to  calculate  the  effect  of  the 

lateral  interactions.  His  assumption  was  that  the  molecules 
t  h 

in  the  i —  layer,  say  n^  of  them,  serve  as  sites  for  the 

-f- 

ni+l  m°lecules  of  the  (i  +  l)s  layer,  and  that  the  n^ 
molecules  are  distributed  randomly  over  the  N  sites  of  the 
lattice  even  though  not  all  N  sites  are  available  to  them. 

The  n^  molecules  of  the  next  layer  are  distributed  randomly 
over  the  n^  sites  available  to  them,  and  hence  are  random 
over  all  N  sites  of  the  lattice.  The  number  of  nearest 
neighbour  pairs  in  each  layer  is  then  calculated  on  the 
basis  of  this  randomness  and  used  in  the  partition  function. 

From  the  partition  function  Hill  found  n  equations 
( A  =  number  of  layers)  in  n  unknowns  (the  0^  Ts  -  the  coverages 
of  each  layer).  The  final  isotherm  is  the  sum  of  the  G±' s. 
Numerical  methods  had  to  be  used  to  find  solutions  for  the 
isotherm.  From  his  calculations.  Hill  found  isotherms  which 
he  says  would  indicate  that  for  sufficiently  large  lateral 
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interactions  and  unrestricted  adsorption  (n  -►  °° ) }  steplike 
isotherms  would  result.  Although  he  did  not  make  calculations 
for  the  n  =  oo  case,  his  isotherms  for  n  >  3  and  large  inter¬ 
action  energies  contained  loops  which  when  cut  according  to 
the  rule  of  equal  areas  showed  that  each  layer  filled 
critically.  In  conclusion  Hill  noted  that  these  results 
followed  from  a  very  crude  approximation  and  that,  although 
it  does  not  happen  in  the  monolayer  case,  the  use  of  a 
higher  order  approximation  or  an  exact  method  might  change 
the  qualitative  nature  of  the  isotherms. 

Halsey  and  Champion  (27)  have  extended  HillTs 

treatment  to  include  the  transmission  of  energy  from  the 

t  h 

adsorbent  out  to  the  n  layer.  E^  is  the  energy  of  attraction 

between  the  first  layer  and  the  surface,  and  Halsey  considers 

that  this  is  transmitted  outward,  decaying  with  the  cube 

of  the  distance  from  the  surface,  so  that  E^  =  E^/nJ  is 

t  h 

the  energy  of  attraction  between  the  n  layer  and  the  surface. 
The  effect  of  this  is  to  spread  the  loops  (i.e.  the  steps) 
of  the  isotherm  on  the  horizontal  scale.  All  of  the  steps 
occurred  at  the  saturation  pressure  (p  =  p  )  in  Hillts 
model,  whereas  in  Halsey*s  model  the  first  five  steps  occur 
at  approximately  p  =  0.28  pQ,  p  =  0.67  PQ.>  P  =  0.86  pQ, 

P  =  0.93  pQ,  and  p  =  O.96  PQ.  The  steps  after  the  first 
three  occur  sufficiently  close  together  that  one  would  not 
expect  to  see  them  in  an  experimentally  determined  isotherm, 
even  if  the  system  were  very  nearly  ideal.  However  the  initial 
steps  are  very  distinct. 
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The  best  example  of  an  experimentally  determined 
stepwise  isotherm  is  the  one  reported  by  Singleton  and 
Halsey  (13).  For  the  system  krypton  on  carbon  black  which 
had  been  graphitized  at  2700°  K,  the  isotherm  shows  two 
very  distinct  steps  in  the  region  1  <  6  <  2  and  2  <  G  <  3  . 
It  has  been  shown  (l4)  that  the  adsorbent  used  had  extremely 
large  crystallites,,  and  thus  was  very  close  to  a  uniform 
surface.  The  fact  that  most  experimental  systems  yield 
smooth  rather  than  step-like  isotherms,  together  with  cal¬ 
culations  which  show  that  the  requirements  necessary  to 
remove  the  steps  in  his  theory  are  entirely  unphysical,  is 
enough  to  make  Halsey  assert  that  a  smooth  isotherm  is 
sufficient  to  conclude  that  the  adsorbent  in  a  system  is 
far  from  uniform.  The  step-like  isotherm  is  what  one  would 
expect  in  an  ideal  system. 

For  a  more  thorough  review  of  these  and  other 
theories  of  adsorption,  the  reader  is  referred  to  the  book 
by  Young  and  Crowell  (28). 
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Chapter  3«  The  Ising  Model  and.  Monolayer  Adsorption 

3 . 1  The  Ising  Model 

In  this  chapter  we  will  discuss  the  work  of 
Ballentine  and  Betts  (8,9)  in  which  they  apply  the  Ising  model 
of  ferromagnetism  to  the  problem  of  localized  monolayer 
adsorption . 

When  Ising  (3)  first  proposed  the  model  in  1925  he 
considered  an  ideal  crystal  of  ferromagnetic  material,  and 
with  each  lattice  point  in  the  crystal  he  associated  a  spin 
variable  ot  which  took  on  the  discrete  values  +1  and  -1 
corresponding  to  "up"  and  "down"  spin  respectively.  He  neg¬ 
lected  interactions  between  atoms  in  the  lattice  which  were 
not  nearest  neighbours  and  he  assumed  that  the  n-n  atomic 

interaction  could  be  expressed  as  -J  a. a.  ,  where  i  and  j 

J 

denote  n-n  spin  sites.  Thus  the  interaction  energy  of  parallel 
spin  pairs  is  -J,  and  that  of  antiparallel  pairs  is  +J.  The 
interaction  of  each  spin  with  the  external  magnetic  field  is 
-H  cu  ,  where  a  spin  ct^  =  +1  is  aligned  in  the  direction  of 
the  magnetic  field,  and  a±  =  -1  is  aligned  opposite  to  the 
field.  The  magnetic  moment  per  spin  has  been  set  equal  to 
one  in  the  above  definition.  The  total  energy  of  the 


y 

The  abbreviation  "n-n"  will  sometimes  be  used  for  'nearest 
neighbour. " 
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ferromagnet  is  then 


E  =  -J  2  a.  o.  -  H  2  or.  , 
n-n  1  J  p  1 

pairs 


and  the  partition  function  is 


(3.1) 


z 


exp 


KT 


2 

n-n 


cr.  o .  + 

i  J 


_H 

KT 


2  a.  } 
i  1 


(3.2) 


following  Ballentine's  notation 

sum  over  all  spin  states  2 

a1=±l 

lattice  of  N  spins. 


where 


represents 


the 
for  a 


Ising  himself  was  only  able  to  solve  the  problem 
for  the  one-dimensional  case.  This  system  did  not  exhibit 
any  critical  behaviour  so  attention  was  turned  to  the  two- 
and  three-dimensional  problems.  Many  workers  found  approxi¬ 
mate  solutions  to  the  two-dimensional  problem,  all  of  which 
exhibited  a  finite  jump  in  specific  heat  at  some  non-zero 
critical  temperature  (Curie  point).  The  first  exact  result 
obtained  for  the  two-dimensional  problem  came  in  1941  when 
Kramers  and  Wannier  (29)  showed  that  the  partition  function 
could  be  derived  by  finding  the  largest  eigenvalue  of  a 
matrix.  For  the  simple  quadratic  lattice  they  found  a 
unitary  matrix  transformation  which  established  a  symmetry 
relation  between  high  and  low  temperature  expressions  for 
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the  partition  function.  From  this  they  deduced  the  location 
and  nature  of  the  singularity  in  the  specific  heat  curve. 

An  excellent  review  of  this,  and  indeed  of  all  the  work 
done  on  the  Ising  model  up  to  i960,  is  found  in  a  paper  by 
Domb  (30).  Other  good  reviews  of  the  Ising  model  are  given 
by  Newell  and  Montroll  (31)  and  Hill  (32).  In  1963 
Fisher  (33)  published  a  shorter  but  more  up-to-date  review 
of  the  problem. 

In  1944  Onsager  (34)  succeeded  in  solving  exactly 
the  zero  magnetic  field  case  of  the  two-dimensional  square 
lattice  Ising  model.  To  date  neither  the  non-zero  field  nor 
the  three-dimensional  cases  have  been  solved  exactly.  In 
1949  Kaufman  (35)  simplified  Onsager’s  matrix  solution.  In 
1952  Kac  and  Ward  (36)  gave  a  combinatorial  solution  based 
on  the  counting  of  a  certain  type  of  graph  on  the  appropriate 
lattice.  Their  solution  was  not  completely  justified  until 
i960  when  Sherman  (37)  showed  that  their  solution  was 
equivalent  to  that  of  Onsager.  In  the  meantime  these  various 
solutions  had  been  extended  to  the  triangular  and  hexagonal 
lattice  cases.  Yang  (38)  extended  the  matrix  method  to 
calculate  the  spontaneous  magnetization  of  a  simple  quadratic 
ferromagnetic  system,  and  this  too  was  later  extended  to 
include  the  triangular  and  hexagonal  lattice  cases. 
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The  most  fruitful  approach  to  the  three-dimensional 
and  non-zero  field  two-dimensional  problems*  for  which  no 
closed  solutions  have  yet  been  found*  is  that  of  exact 
series  expansions.  Comparisons  of  the  approximate  methods 
with  exact  expansions  have  shown  that  the  series  method  can 
always  be  made  more  accurate  than  any  approximation  without 
having  to  calculate  an  excessive  number  of  terms.  This  and 
the  fact  that  the  approximations  distort  the  results  in  the 
critical  region*  make  the  exact  series  approach  a  more 
profitable  avenue  to  pursue*  at  least  for  the  present. 

Van  der  Waerden  (39)  and  later  Wakefield  (40) 
used  the  combinatorial  technique  of  graph  counting  to  obtain 
the  high  temperature  series  expansion  of  the  partition 
function  for  zero  field.  Their  expansion  coefficients 
involved  the  counting  of  closed  graphs.  Oguchi  (4l)  derived 
a  series  for  the  zero-field  magnetic  susceptibility  on  a 
cubic  lattice*  which  involved  counting  all  graphs  with  two 
odd  vertices.  On  the  other  hand  Domb  (42)  used  a  matrix 
method  to  derive  the  series  for  the  partition  function  on 
a  simple  quadratic  lattice.  Recently  Domb  and  Sykes  have 
done  extensive  work  in  extending  the  graph  counting  techniques. 
These  techniques  are  described  in  detail  in  chapter  ol 
the  paper  by  Domb  (30)*  and  were  used  in  calculating  the 
series  expansions  which  occur  in  Chapter  4  ol  this  thesis. 
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3.2  Application  to  Localized  Monolayer  Adsorption  with 

Lateral  Interactions 

In  applying  the  Ising  model  to  adsorption  we 
consider  the  adsorbent  to  be  a  regular  two-dimensional  lattice, 
the  vertices  of  which  are  all  possible  sites  on  which 
molecules  of  the  adsorbate  may  settle.  Thus  each  site  may 
be  in  one  of  two  states  denoted  as  follows:  cj  =  +1  if 
the  site  is  occupied  by  an  adsorbed  molecule,  and  o\  =  -1 
if  the  site  is  unoccupied.  cu  is  now  an  "occupation" 
variable,  although  from  time  to  time  we  may  refer  to  it  as 
"spin"  in  analogy  to  the  ferromagnetic  case.  The  subscript  i 
simply  refers  to  the  lattice  site  and  may  take  the  values 

2,  .  .  .  ,  N.  We  can  write  the  total  number  of  adsorbed 

molecules  as 

N  =2  (1  +  o. )/2  (3.3) 

a  .  i 

i 

and  the  total  number  of  n-n  pairs  of  adsorbed  molecules  as 

N  =  2  (1  +  (T.  )(l  +  a  )/4  (3-4) 

p  n-n  1  J 

The  model  includes  an  attraction  between  every 
adsorbed  molecule  and  the  surface,  characterized  by  an 
interaction  energy  -E  ,  and  an  attraction  between  adsorbed 
characterized  by  an  interaction  energy  -V. 


n-n  pairs. 
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Including  the  chemical  potential  \i  of  the  system  the  energy 
may  be  written  as 


E  =  -(u.  +  Es)  Na  -  V  Np  .  (3.5) 

Using  (3. 3)  >  ( 3  -  *4 ) ,  and  (3-5)  the  grand  partition  function 
of  the  adsorbed  system  is 


.  \i  +  E 

2  exP  ( — p — “  2  ( 1  +  a  )+  j-  2  (l+cr  )  (l+a  ) } .  (3*6) 

( cr.  )  i  n-n  J 

v  1 

Expanding  the  product  in  the  second  term  of  (3.6)  and 
remembering  that  there  are  n-n  pairs  so  that 


2  a.  =  £  2  o.  , 
1  2  .  l 

n-n  i 


(3.7) 


we  get 


7 


=  exp  +  cr-  +  ^))  2  exp  ( 

(^) 


KT  v2  '  2  '  8 


(E  +M.+SI) 

f  l  Vw  2  Vj  ] 

2  KT  i  n-n 


If  we  let  H  =  i(E  +  d  +  7p)  and  J  =  ^  ,  then 


7,= exp  [k!  (2  +  “I  +  ]  z> 


(3.8) 


where  Z  is  the  standard  Ising  partition  function  of 
equation  (3.2).  The  coefficient  of  Z  in  (3*8)  is  a 
function  of  temperature.,  but  otherwise  it  is  independent  ol 
the  state  of  the  adsorbed  system. 
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3 • 3  Series  Expansions  of  the  Partition  Function 

We  wish  to  have  expansions  of  Z  in  both  high  and 
low  temperature  variables  -  i.e.  variables  which  tend  to 
zero  as  the  temperature  tends  to  infinity  and  zero  respectively. 
For  a  given  configuration  of  the  system  we  say  that  of  the 
N  spins,  n  of  them  are  "down"  or  negative  spins,  and  of 
the  ^  n-n  pairs,  m  of  them  are  antiparallel.  Then  we  can 
write 


and 


Z  cr.  =  N  -  2n  , 

l 

l 


Z  cr.  cr .  = 
i  J 
n-n  0 


qN 

2 


-  2m 


Bringing  the  factors  which  are  independent  of  the  configura¬ 
tion  of  the  system  to  the  left  of  Z  we  have 


Z 


rNH  ,  qNJ 
exP  '-jqp  2KT 


exp 


f  2nH 
1  KT 


2mJ  •) 
KT  J 


We  now  choose  the  variables 

-2J/KT 

x  =  e  and 


-2H/KT 

y  =  e 


(3.9) 


noting  that  x  —  0  and  y  —  0  as  T  —  0.  In  the  sum  over 
states  Z  ,  the  number  of  negative  spins  n  will  take 
on  all  ^  ^  values  from  0  to  N,  and  the  number  of  anti- 
parallel  n-n  spin  pairs  m  will  take  on  all  values  from 
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0  to  qn  for  a  given  n.  Hence  we  can  write 
QN  N  n  qn 

Z=x  2y2Z  2  g(n,  m;  N)  xm  yn  ,  (3.10) 

n=0  m=0 


where  g(n,  m;N)  is  the  number  of  configurations  on  the 
lattice  of  n  negative  spins  such  that  they  form  m  anti¬ 
parallel  spin  pairs.  Equation  (3.10)  is  the  low  temperature 
(L.T.)  expansion  for  Z,  and  the  calculation  of  the  actual 
series  involves  the  counting  of  the  configurations  just 
described. 

The  high  temperature  (H.T.)  series  can  be  found 
by  writing  (3.2)  in  the  form 


Z 


TT 

n-n 


(J/KT) a 
e 


(H/KT) 


a, 


(3.11) 


By  expanding  the  exponential  factors  above,,  using  the 
following  relation 


(a.r  =  (^j)  =  1 

and  grouping  the  terms  of  the  series  properly,  we  can  write 


Tr-q-,0”  .  (7  .  J 

i  1  J  =  (l  +  ct.  a.  u)  cosh  KT 

1  J 


and 


H 


(3.12) 


a 


KT  1  =  (l  +  a.  v)  cosh 


H 

KT 
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where  u  _  tanh  ^  and  v  =  tanh -[Jp  .  Substituting  the 
relations  (3-12)  into  equation  ( 3 . 11 )  gives 


qN 


J 


z=  (cosh  *  (cosh  §=■)  2  TT  (l+o.  o.  u)  JJ  (l+d,  v).  (3.13) 


N 


KT 


KT 


(cn)  n-n 


i  J 


i 


We  expand  the  products  in  (3*13)  and  consider  the  coefficient 
m  £ 

of  v  u  .  It  will  be  the  sum  of  all  possible  products  of 

i  different  factors  of  the  form  a. a.  and  m  different 

i  J 

factors  of  the  form  cl .  Each  such  product  may  be  represented 

graphically  if  we  let  a .  cr .  be  represented  by  a  line  or 

a  J 

t  h  t  h 

bond  between  the  i —  and  j —  sites  on  the  lattice.,  and 

t  h 

be  represented  by  a  circle  around  the  k —  site.  The 

m  H/ 

coefficient  of  v  u  then  is  represented  by  the  sum  of  all 
the  graphs  of  l  lines  and  m  circles.  However  in  the  sum 
over  states,,  if  a  particular  term  contains  any  factor  cl 
to  an  odd  power,  its  contribution  will  be  zero.  Hence  graphs 
which  contain  only  vertices  which  either  are  uncircled  and 
are  the  meeting  point  of  an  even  number  of  lines,  or  are 
circled  and  are  the  meeting  point  of  an  odd  number  of  lines, 
will  be  the  only  ones  to  contribute  to  the  final  expansion. 
The  form  of  the  H.T.  expansion  for  Z  is 


SE  H  N 

Z  =  (cosh  )  2  (cosh  £tjt) 


2X 


w,  z  2  Jj{t,  m;N)  vV 

^  jg=0  m=0 


y 


(3.14) 
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where  (I,  m;  N)  is  the  number  of  graphs  of  Jb  lines  which 
contain  m  odd  vertices  and  can  be  drawn  on  a  lattice  of 
N  sites.  For  large  values  of  Jb  the  upper  limit  of  2 Jb 
on  m  will  never  be  realized,,  but  for  small  Jb  a  graph 
of  Jb  lines  will  have  at  most  2 jb  vertices. 


3.4  Calculations  Using  the  Partition  Function 


Since  we  are  considering  a  regular  lattice,  in 
which  all  sites  must  be  identical,  we  will  let  the  number 
of  sites  N  become  infinite,  or  consider  the  lattice 
boundaries  to  be  cyclically  connected.  In  this  way  we 
avoid  any  edge  effects  which  might  otherwise  have  entered. 


However,  we  must  first  see  how  the  partition 
function  may  be  used  to  derive  the  physical  quantities  we 
are  interested  in.  From  statistical  mechanics  we  know  that 
we  can  find  equilibrium  values  of  a  physical  quantity  by 
averaging  that  quantity  over  all  the  possible  states  of 
the  system.  For  instance  the  internal  energy  of  the 
ferromagnetic  system  has  an  equilibrium  value  given  by 


2  E  e 
s 

E  =  - - 


-Eg/KT 


2  e 
s 


-ES/KT  ' 


(3.15) 


where  E  refers  to  the  energy  of  a  particular  state  of  the 
s 

system,  and  2  is  the  sum  over  states. 

s 


From  this  it  is 
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easy  to  show  that  the  internal  energy  per  spin  is 


E  =  KT"  It  (!  ln  z)  • 


(3.i6) 


The  specific  heat  is  simply  obtained  by  taking  the 

derivative  of  E  with  respect  to  T.  The  magnetization  per 
spin  can  be  found  in  a  way  analogous  to  (3.15)  to  be 


1 


M  =  KT  3h  %  ln  z)  • 


(3.17) 


Ballentine  has  shown  that  the  fractional  coverage  9  of  the 
surface  in  the  adsorption  problem  is  related  to  M  by 


0  =  |  (M  +  1)  . 


(3.18) 


If  we  remember  how  x,  y,  u,  and  v  were  defined  in  terms  of 
J,  H,  and  T  we  can  write  and  9  as  functions  of  these 

variables . 

Physically  we  know  that  lim  ^  ln  Z  must  be 
finite.  Thus  the  series  for  ln  Z  must  be  independent 

of  any  powers  of  N  greater  than  the  first.  Ballentine  has 
rigorously  derived  the  ln  Z  series  from  the  series  for  Z, 
and  used  this  condition  on  the  coefficients  of  the  ln  Z 
series  as  a  check  of  his  work.  It  can  be  shown  that  the 
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coefficient  of  a  given  term  in  the  In  Z  series  is  simply 
the  term  linear  in  N  in  the  polynomial  coefficient  of  the 
corresponding  term  in  the  Z  series.  For  example,  in  the 
H.T.  Z  series,  the  coefficients  m;  N)  are  polynomials 

in  N,  and  the  coefficient  of  the  vmu^  term  in  the  In  Z 
series  is  the  term  linear  in  N  in  h  U,  n)  . 

Using  the  series  for  In  Z,  Ballentine  has 
calculated  the  adsorption  isotherms  given  by  this  model. 

For  temperatures  greater  than  a  critical  temperature  T  , 
the  isotherms  are  smooth  curves  such  that  Q  increases  as 
the  pressure  increases.  However  for  temperatures  below  Tc 
there  is  a  finite  discontinuity  in  the  isotherm,  occurring 
at  a  pressure  P  .  Using  the  relation  between  the  chemical 
potential  \i  and  the  pressure 


U  2  3 

KT  _  /  h  \  2  P_ 
"  12t r  mKT  '  KT 


(3.19) 


for  an  ideal  gas  (when  the  gas  phase  is  in  equilibrium  with 
the  adsorbed  phase,  the  chemical  potentials  of  the  two 
phases  are  equal),  one  can  show  that  this  critical  pressure  Pq 
corresponds  in  this  model  to  the  case  H  =  0.  Yang  (38) 
has  calculated  an  exact  formula  for  the  spontaneous  magneti¬ 
zation  (M  when  H  =  0)  of  the  simple  quadratic  Ising  ferro- 
magnet.  Ballentine  has  used  this  formula  to  calculate  exact 


values  for  the  isotherms  right  at  the  discontinuity.  The 
critical  temperature  for  this  system  on  the  simple  quadratic 
lattice,  first  found  by  Kramers  and  Wannier  (29)  is  given  by 

-  2J 
KT 

X  =e  c  =  v/2  -  1  =  0.4142  .  (3.20) 

3. 5  Two  or  More  Adsorbed  Layers 

Ballentine  has  also  considered  an  extension  of 
this  model  to  treat  two  or  more  layers  of  adsorbed  molecules. 
However  this  model  is  based  on  the  physically  unrealistic 
premise  that  the  state  of  a  site  in  a  higher  layer  is 
independent  of  the  states  of  the  sites  directly  below  it. 

This  allows  gas  molecules  to  be  adsorbed  on  top  of  holes  in 
a  lower  layer.  However,  Ballentine  studied  this  model  because 
it  was  a  simple  picture  of  cooperative  phenomena  and  still 
displayed  a  change  in  critical  point  with  an  increasing 
number  of  layers.  Ballentine  has  shown  that  the  treatment 
of  the  n-layer  case  for  2  <  n  <  «>  is  very  difficult.  For 
n  =  1,  2,  oo  layers  of  a  simple  cubic  lattice  structure  we 
are  dealing  with  homogeneous  lattices,  each  site  having 
4,  5,  and  6  nearest  neighbours  respectively.  However  for  any 
other  value  of  n  the  lattices  are  inhomogeneous,  sites  in 
the  top  and  bottom  layers  having  five  nearest  neighbours 
while  interior  sites  have  six.  Ballentine Ts  method  for 
calculating  the  series  coefficients  is  not  as  successful  for 
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inhomogeneous  lattices  as  it  is  for  homogeneous  ones.  He 
has,  however,  made  an  accurate  calculation  of  the  zero- 
field  specific  heat  together  with  a  good  estimate  of  the 
critical  temperature  for  the  two-layer  case.  For  the  three 
and  four  layer  cases  he  has  made  good  calculations  of  the 
specific  heat  away  from  the  critical  region  and  has  made 
rough  estimates  of  the  location  of  the  critical  point.  For 
the  limiting  case  of  an  infinite  number  of  layers,  the  zero- 
field  specific  heat  and  the  critical  temperature  can  be 
found  from  WakefieldTs  (40,43)  exact  series  expansion  for 
the  simple  cubic  lattice.  Ballentine’s  results  show  that  there 
is  a  smooth  variation  in  critical  temperature  with  the  number 
of  layers.  This  is  shown  in  fig.  1,  which  is  reproduced  from 
BallentineTs  thesis. 


Fig.  1 


The  variation  in  critical  point  with 
the  number  of  layers  for  the  Ising  model. 
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Chapter  4.  The  Model  for  Two-Layer  Adsorption  with  Lateral 

Interactions 


4 . 1  The  Model 

We  saw  in  the  last  chapter  that  the  two-layer  Ising 
model  when  applied  to  adsorption  permits  the  possibility  of 
a  second  layer  molecule  being  adsorbed  on  top  of  a  hole  or 
vacancy  in  the  first  layer.  One  of  the  principal  reasons  for 
suggesting  a  new  model  is  to  eliminate  this  unrealistic 
possibility.  The  model  which  we  propose  is  as  follows.  We 
consider  the  adsorbent  to  be  a  two-dimensional  square  lattice 
of  adsorption  sites.  Each  site  is  capable  of  containing 
zero,  one,  or  two  adsorbed  gas  molecules.  The  first  molecule 
on  a  given  site  is  said  to  be  in  the  first  layer,  and  the 
second  one  in  the  second  layer.  To  each  site  we  assign  an 
occupation  ("spin")  variable  cn  which  takes  on  the  values 
-1,  0,  and  +1  corresponding  to  occupancies  of  0,  1,  and  2 
molecules  respectively.  We  say  that  there  is  an  interaction 
energy  -E  between  molecules  in  the  first  layer  and  the 
surface  of  the  adsorbent,  an  interaction  energy  -JH  between 
nearest  neighbour  pairs  of  molecules  in  the  same  layer,  smci 
an  interaction  energy  -Jy  between  two  molecules  adsorbed 
on  the  same  site  (and,  by  assumption,  in  difierent  lasers). 

We  assume  that  all  interactions  are  sufficiently  short  range 
that  we  may  neglect  all  but  n-n  interactions. 
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This  model  should  not  be  confused  with  the  spin  1 
Ising  model  for  a  ferromagnet  composed  of  atoms  with  three 
atomic  spin  states.  The  interaction  between  molecules  on 
the  same  site  in  the  adsorption  problem  has  no  analogue  in 
the  ferromagnetic  problem.  Thus  there  will  not  be  the  same 
correspondence  between  the  partition  functions  for  the  two 
cases  that  there  was  in  BallentineTs  work. 


We  can  write  the  following  expressions  using  the 


occupation  variable  cu  : 


Number  of  molecules  in  first  layer 


Number  of  molecules  in  second  layer 


N  cr  (l  -  ct  ) 

=  2  (l  +  -i__ - —  )  , 

i=l  ^ 

Mi  +  M  % 

=  2  (— —  s - —  )  - 

i 


Total  number  of  molecules  =  2  (l  4-  cr  ) 

i 


(4.1) 


Using  the  relations  (4.1)  we  can  write  down  the  total  energy 
of  the  adsorbed  system. 


E  =  -|jl  [2  (1  +  o±)] 

cr,  ( 1  -  a  ) 


i 


-E  [2  (1  + 

s  . 

l 


Ml 


o.fi-ffj  cr.(i-a.)  a  (1+C7  )  cr  (l+a  ) 

-j  [z  ((l  +  1  )(1  +  J  J--)  +  ( - 2 - } 

n-n  ^  d 

a.  (l  +  cr.  ) 

-JV  - 2 - • 


2 


(4.2) 


We  expand  this  expression  and  use  equation  (3-7)  to  get  the 
partition  function  Z. 
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Z  = 


(a1=0,±l) 


pxd  ftU  + 
p  lKT  KT 


E  N  qJ„N 


+ 


H 


2KT 


) 


2KT  n^p_  (cJi'7j  +  Ji2(Jj2' 


(4.3) 


E 


J, 


,  ,  _s_  3V_ 

lKT  2KT  2KT 


qJu  -E  jT7 

)  2  °4  +  (oW  +  c 


2KT '  7  ~i 

l 


2KT  2KT 


Q  J  u  o 

2kt)  2  cr.  } 


There  are  three  parameters  in  this  theory  which  are 

completely  arbitrary;  E  ,  J^,  and  JTT.  However  because  (4.3) 

is  sc  cumbersome  and  an  expansion  of  it  as  it  stands  would 

be  very  difficult,  we  will  simplify  (4.3)  by  requiring  that 

2 

the  coefficients  of  the  2  a.  and  2  a.  terms  be  zero. 

•  1  1 
i 

JV  "  Es  ”  q  JH  “  0  * 

JTr  +  E  +  q  J„  +  2\j,  0  .  (4.4) 

The  first  equation  in  (4.4)  essentially  restricts  the 
physical  systems  to  which  we  may  apply  the  model  as  it  now 
stands.  The  two  equations  (4.4)  reduce  the  number  of  free 
interaction  parameters  from  three  to  one  and  any  system  on 
which  this  is  used  must  be  such  that  the  interactions  are 
in  the  ratio  which  (4.4)  requires.  The  presence  of  q. 
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in  the  second  equation  means  (see  equation  (3.19))  that  we 
have  a  fixed  p  -  T  condition  for  this  special  case,,  and 
that  we  will  be  unable  to  calculate  adsorption  isotherms 
as  long  as  (4.4)  holds.  The  partition  function,  then,  in 
which  we  are  interested  is 

z  =  exp  [Jjjr  (n  +  E  +  -Lii)]  2  exp[  Uf  2  (5,0,+  ^)] 

^  ad  (ai=0,±l)  2KT  n-n  1  J  1  J  • 

(4.5) 


4 . 2  The  low  Temperature  Series  Expansion  for  Z 

For  the  purpose  of  finding  a  L.T.  expansion  we 
write  (4.5)  as 


z  =  [  A  ( T )  ] 


N 


(ffi 


2 

:0, 


.  TT 

±1 )  n-n 


exp 


hi 

1 2KT 


^ i°' j  + 


a.  2a.2)  ] 

1  J  1  9 


(4.6) 


We  can  see  that  the  part  of  the  energy  of  the  adsorbed  system 
which  is  dependent  upon  the  configuration  of  the  systemic,  is 


J 


P  =  2  [  (cr.  cr.  +  o^a*)  ] 

C'  L  '  i  .1  i  .1  ' 


2  2 


n-n 


2  i  J  i  J 


(4.7) 


2  2 


By  examining  the  possible  values  of  the  function  uc.  +  ■  . 


below 
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\i 

-1 

0 

+1 

a . 

-1 

+2 

0 

0 

0 

0 

0 

0 

+1 

0 

0 

+2 

2  2 

f  (cr.,a. )  =  a.  a.  +  a.  a. 
i’  J  i  J  i  j 


we  can  see  that  £  will  have  its  minimum  value  when 

either  all  "spins"  are  +1  or  all  are  -1.  The  ground  state 
is  doubly  degenerate  and  its  energy  is 


For  a  higher  energy  state  the  configurational  energy  is 

incremented  by  for  every  n-n  pair  of  spins  of  type  0 

n 

(the  pair  (a., a.)  is  type  0  or  2  depending  upon  whether 

f (a  ,  o.)  =  0  or  2)  which  is  introduced  into  the  configuration 
v  i  J 


The  sum  over  states  2  includes  every 

(ct  =0,±1) 

possible  configuration  of  the  system,  so  the  L.T.  expansion 

will  start  with  the  ground  state,  and  each  succeeding  term 

will  represent  the  successively  higher  energy  states 

corresponding  to  the  formation  of  more  and  more  spin  pairs  of 

-Jh/KT 

type  0.  We  choose  the  variable  x  =  e  ,  and  write 

the  L.T.  series  as 


V 


43. 


z 


[a(t)]N  exp  [-7^] 


qNJ 


H' 


2KT 


qN 

2 

2 

n=0 


2g(n) 


x 


n 


(4.8) 


2g(n)  is  equal  to  the  number  of  configurations  on  the  lattice 
which  contain  n  "spin"  pairs  of  type  0.  We  have  chosen 
to  write  this  as  2g(n)  because  every  configuration  in 
which  we  are  interested  is  related  to  one  other  similar 
configuration  with  the  same  energy  which  can  be  obtained 
by  the  transformation  ^  .  Thus  in  counting  the 

configurations  to  evaluate  g(n)  we  agree  to  count  only 
configurations  which  have  a  prodominance  of  +1  "spins". 

To  aid  in  the  counting  of  the  configurations  g(n) 
we  consider  the  following.  In  a  given  configuration  of  spins 
on  the  lattice  we  draw  a  line  to  separate  each  "spin" 
pair  of  type  0.  These  lines  will  form  a  "graph"  on  the 
reciprocal  lattice.  Since  the  lattice  is  cyclically  connected 
all  graphs  will  be  closed.  Hence  to  calculate  g(n)  we  must 
count  all  closed  graphs  of  n  bonds  on  the  reciprocal 
lattice,  the  weight  or  contribution  of  each  graph  being  the 
number  of  ways  we  can  arrange  +1,  0,  and  -1  spins  in  the 
interior  of  each  lattice  cell  (the  number  of  +1  "spins"  must 
be  greater  than  the  number  of  -1  spins  )  so  that  pairs  oj 
type  0  are  always  separated  by  a  line  and  pairs  of  type  2  are 
not.  Since  the  simple  quadratic  la^^ice  is  self -reciprocal 
we  will  be  interested  in  closed  graphs  on  a  simple  quadratic 


lattice . 
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],□) 


Example :  The  separated  closed  configuration  ( 

of  13  bonds  can  be  drawn  on  the  s.q.  lattice  in  4N2  -  60N 
ways.  Each  of  these  ways  makes  a  contribution  of  2, 
corresponding  to  the  following  ways  of  "filling”  the  graph: 


+  +  +  +  +  +  + 
3  +  I°l  + 


+ 


+  +  +  +  +  +  + 


+  +  +  +  +  +  + 
+ 


o]+  3  + 


+++++++ 


Thus 
2. 


( 


-□>  makes  a  total  contribution  of 


8  N  -  120  N  to  g (13 )  . 


4 . 3  The  High  Temperature  Series  Expansion  for  Z 

To  make  the  H.T.  series  expansion  we  consider  Z 

in  the  form  (4.6)  and  expand  the  exponential  term  in 
2  2 

a.  cr.  +  a.  a.  in  an  infinite  series.  We  let 
1  J  1  J 


Pij  ~  CTiaj  +  ai  °*j  ’ 


(4.9) 


and  we  can  easily  prove  that  for  the  values  01  =  0,  ±1 


p.  .  =  2  p.  . 

iJ 

n  __  p(n-l) 

pij  =  2  pu 


(4.10) 


The  expansion  becomes 

z  =  [a(t))n  2  x  T7  (1+St  pij+  li  4fr)2  2Pij+l!^T)322Pij+ 

(a  =0,±l)  n-n 


•  0  • 
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Rearranging  the  terms  we  get 


z  =  [a(t)]n  2  tt  U+^ii[i+5i  +1  (hi)2-!.!  (ti) 3  , 

(a.=0,±l)  n-n  2  u+kt+?.W  'p.  W  + 


2  L  KT  2.'Ut'  '  3J  'KT 


]- 


P 


And  finally 


N 


z  =  [A(T)r  2 


e5i  -  1 

KT  J- 


(a1=0,±l) 


^  +  Pif'  2 


)  )  •  (4.11) 


n-n 


We  choose  the  H.T.  variable  to  be 


w 


e5i  -  1 
=  = 


1  ,  -1 
O  (X 


-  1) 


and  note  that  as  T  ,  w  0. 

When  we  expand  the  product  in  (4.1l)  we  note  that 
n 

the  coefficient  of  w  will  be  the  sum  of  all  possible  products 

T 

of  n  p.  .  s  .  We  now  associate  with  each  p.  .  a  bond  on 
ij  ij 

the  lattice  drawn  between  the  n-n  sites  i  and  j.  Let  us 

now  consider  the  contribution  that  a  given  configuration  of 

n 

n  bonds  will  make  to  the  coefficient  of  w  when  we  sum  over 

all  states  of  the  system.  The  sum  2  will  contribute 

ak  «0,±1 

a  factor  of  3  if  the  subscript  k  is  absent  from  all  the 
p  ?s  in  the  product  (i.e.,  if  no  line  goes  through  the  k— 
vertex  of  the  lattice).  Thus  if  m  vertices  have  lines 
passing  through  them,  the  N-m  "empty"  vertices  will  make 
a  contribution  of  e^N-m^  in  the  sum  over  states.  The 


■H  OJ 


46. 


contribution  from  the  m  vertices  which  are  touched  by 
lines  on  the  graph  will  be  in  two  parts.  A  factor  of  2 
enters  for  each  such  vertex,,  and  in  addition  a  factor  of 
2  enters  for  each  "irreducible  loop"  formed  by  the  graph. 

By  "loop"  we  mean  a  closed  path,  through  all  or  part  of  a 
graph,  around  which  one  may  travel  and  return  to  the  starting 
point  without  retracing  a  step.  By  "irreducible  loop"  we 
mean  the  smallest  loop  which  returns  to  a  given  starting 
point.  An  "irreducible  loop"  may  have  a  boundary  in 
common  with  another  "irreducible  loop"  but  it  may  not 
circumscribe  any  part  of  the  area  circumscribed  by  another 

contains  the 


one.  For  example  the  configuration 


three  irreducible  loops  shown  in  fig.  2(a).  The  loops 
shown  in  figs.  2(b)  -  2(d)  for  the  same  configuration  are 
not  all  irreducible. 


' - 

l 

r  ■  ■ 

3 

2 

1 


\ 

— 

2I 

- 

— S 

(a)  (b)  (c)  (d) 

Pig.  2.  (a)  The  three  "irreducible  loops"  for  the  graph 

(b)  &  (c)  -  loop  1  is  not  a  minimal  path. 

(d)  loop  1  circumscribes  the  whole  area  of  loop  2 


If  there  are  p  irreducible  loops  in  a  given  graph  of  n  bonds, 
its  total  contribution  to  the  coefficient  of  wn  in  the  sum  over 
states  will  be  3(N_m)2 .  The  following  relation  for  the  pL .  s 
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2 

ar°'±1 


(4.12) 


together  with  (4.10)  makes  it  easy  to  show  that  the  above 
analysis  is  correct. 

If  we  bring  a  factor  3^  outside  the  summations 
we  can  write  the  H.T.  expansion  for  Z  as 

N 

Qg  2n 

z  =  [3A(T)]N  2  2  2  (|)ra  2p  G(n,m,(p);  N)  wn,  (4.13) 

n=0  m  p  ^ 

where  G(n*m,(p);  N)  is  the  number  of  graphs  of  n  bonds 
which  touch  m  vertices  and  form  p  irreducible  loops 
which  can  be  drawn  on  a  lattice  of  N  sites.  The  counting 
of  the  graphs  required  for  G(n*m* (p) ;N)  is  quite  laborious 
because  all  possible  configurations  of  n  bonds  are  involved 
and  each  configuration  must  further  be  classified  according 
to  the  number  of  vertices  it  touches  and  the  number  of 
irreducible  loops  it  includes  before  its  contribution  can 
be  calculated.  However*  as  we  shall  see  in  the  next  section 
these  configurations  have  been  counted  and  classified  for 
small  values  of  n. 


4.4  The  Calculation  of  the  Various  Series 

As  we  have  seen  the  calculation  of  the  series 
coefficients  for  the  two  expansions  involves  the  counting 
of  graphs.  This  method  has  been  used  extensively  by  many 
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people  working  in  the  field  of  cooperative  phenomena. 

Chapter  5  of  the  paper  by  Domb  (30),  written  in  collaboration 
with  Sykes,  gives  a  thorough  treatment  of  the  methods 
available  to  count  virtually  any  type  of  graph.  Most  of  the 
constants  used  to  calculate  the  series  coefficients  were 
taken  from  the  tables  in  Appendix  3  of  Domb’s  paper. 

However  the  closed  graphs  of  order  greater  than  ten,  i.e. 
of  more  than  ten  bonds,  used  in  the  L.T.  series,  and  the 
separated  graphs  of  order  six  used  in  the  H.T,  series  were 
calculated  particularly  for  this  work.  The  equations  and 
checks  used  in  the  calculation  of  the  sixth  order  separated 
lattice  constants  are  given  the  Appendix  of  this  thesis. 

Table  1  shows  how  the  L.T,  coefficients  were 
calculated.  Each  contributing  graph  is  listed  together 
with  the  lattice  constant  for  that  configuration  and  the 
number  of  ways  it  can  be  ’’filled"  according  to  the  prescrip¬ 
tion  in  section  4.2.  Table  2  provides  a  similar  breakdown 
of  the  calculation  of  coefficients  for  the  H.T.  series 


for  Z. 


Table  1:  Breakdown  of  the  calculation  of  the  low  temperature 

series  coefficients  g(n). 


n 

No.  of 
bonds 

Configuration 

Value 

latt. 

const. 

Ways 

to 

fill 

Configuration 

Value 

latt. 

const. 

Ways 

to 

fill 

Coeff . 
of  N  in 

... 

4 

□  ;  p4 

N 

2 

+•2 

6 

O  ’  p6 

2N 

1 

4*2 

7 

O  ;  P7a 

2N 

3 

4*6 

8 

O  5  p8 

7N 

1 

1  1  . — ,P8h  and 

O  P4’P4 

n2-5n 

2 

4 

-3 

9 

c 

^  f  p9k 

12N 

1 

♦12 

10 

O1  Pi° 

28N 

1 

Q 

4n 

5 

( 

) 

2N 

5 

OO  „10. 

8n 

2 

c 

3 

4N 

l 

Q 

p6’p4 

□ 

2N2 

-24N 

2 

+30 

11 

o 

48n 

— — — — 

l 

6 

4n 

1 

<SX> 

8n 

6 

P) 

P?a>P4 

□ 

2N2 

-24N 

6 

-44 

12 

O' 

124n 

1 

0 

0 

12N 

2 

cC- 

6N 

1 

03 

16n 

2 

12N 

1 

8n 

1 

o 

32N 

1 

€ 

3  ■ 

N 

7  ‘ 

-  u  * .. 


Table  1  (continued). 
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Table  2:  Breakdown  of  the  calculation  of  the  high  temperature 

series  coefficients  G(n,m, (p) ;N). 


n 

G(n,m, (p) ;N) 

Lattice  constants  which 
contribute  to  G(n,m,(p);N) 

Coeff.  of  N  in 
G(n,m, (p) ;N) 

1 

G(1 , 2, (0) ;N) 

</) 

42 

2 

G(2,3,(0);N) 

(A) 

+6 

G(2,4,(0);N) 

</,/) 

-7 

3 

G(3>,(0)iN) 

(A/),  (A) 

4  22 

G(3,5,(0);N) 

(A,/) 

-60 

0(3,6, (0);N) 

</,/>/) 

.116 

*  3 

G(4,4, (1 ) ; N) 

(O  ) 

4l 

G(4,5,(0);N) 

(X),  (AA),  (AZ.) 

487 

G(4,6,(0);N) 

(A/,/),  (A,/),  (A  ,A  ) 

-4o6 

G(4,7, (0) ;N) 

(A,/,/) 

+578 

G(4, 8, (0) ; N) 

(/,/,/,/) 

.521 

2 

5 

G(5,5,(l);N) 

(uf ) 

48 

G( 5,6, (1 ) ;N) 

(□,/) 

-12 

G(5,6,(0);N) 

(AA/),  (W_),  (v5), 

<X>>  <X^> 

+364 

G( 5, 7, (0) ; N) 

1 

I 

(/V,/  ),  (A,/) 

-2500 

Table  2  (continued)* 
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n 

G(n,m, (p) ; N) 

Lattice  constants  which 
contribute  to  G(n,m,(p);N) 

Coeff.  of  N  in! 
G(n,ra,(p);N)  j 

G(5,8,(0);N) 

(/V,/,/),(/N, /,/),( A,A  ,/) 

*5822 

G(5,9,(0);N) 

(A ,/,/,/) 

-5644 

G(5,10,(0);N) 

(/,/,/,/,/) 

5 

6 

G(6,6,(1);N) 

+5>+ 

G(6,7,(1);N) 

(Cf  ,/),(□, A) 

-168 

G(6,7,(0);N) 

(AAA ) ,  ( AAZ-  ) ,  (/xK/)  >  (X^ ) 

)  >  <yv> » ()"< ) 

*1574 

G(6,8,(1);N) 

(□,/,/) 

+126 

G(6,8,(0);N) 

(AAA/MWZ-./MX'  ,/) 

<v4  >/mX>/maa,/\) 

(AZ.  ,A  ),(X»A  )><A/>A/) 

<AAA),(A,A> 

-14681 

G(6,9,(0);N) 

(AA,/,/),(AZ,/,/),(X./,/) 

(A/, A,/), (A, A  ,/),(A  ,A  ,A  ) 

>49034 

G(6,10,(0);N) 

(A/, A/,/), (A,/,/,/) 

(A  >A  »/*/) 

-76098 

G(6,11j<0);N) 

(A  >/ >/ >/ >/) 

+560^0 

G(6,12,(0);N) 

47644 

3 
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As  we  said  in  Chapter  3  when  discussing  the  monolayer 

case.,  since  lim  In  Z)  exists  physically  we  know  that 

N  -*■  oo 

all  powers  of  N  greater  than  the  first  must  cancel 

out  when  we  convert  from  the  Z  series  to  the  In  Z  series. 

This  was  checked  for  the  two-layer  case,  using  the  formulas 
for  the  series  In  [f(x)]  in  terms  of  the  series  f(x) 
which  are  given  in  Appendix  I  of  Ballentine's  thesis (8), 
and  found  to  be  correct.  The  coefficient  of  a  power  of  x 
in  the  ^  In  Z  series  is  simply  the  coefficient  of  N  in  the 
coefficient  of  the  corresponding  power  of  x  in  the  Z  series. 

We  will  now  write  down  the  L.T.  and  H.T.  series  for  ~  In  Z. 

p+E  +q  J  00  -I 

L.T.  ^  In  Z  =  - — -  +  2  [coeff.  of  N  in  g(n)]x  +  ^  In  2  . 

^  n=l 

(4.14) 

Because  of  the  conditions  (4.4)  the  first  term  is  zero,  and 
in  the  limit  N  00  the  last  term  becomes  zero.  Hence  we 
have  (using  the  results  of  Table  l) 


lim  b  In  Z 
N 


2xV2x^+6xr',-3x8+12x9+30x10-44x11+2^x12+44x13+  ...  . 

(4.15) 


qj 


H.  T.  ^ln  Z 


=  In  3" 


p+E 


s  2 


H  00 
Z 


KT 


+ 


n=l 


Z  Z(|-)m2p[coeff .of  N  in 
m  p  p 


G(n,m, (p) ;  N) 3 


n 


w 


,  (4.16) 
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and 

lim  i  In  Z 

N  - 


In  3- 


15h  +  ^  32  w2 

KT  9^  BT 


352 

2187 


832 

5551 


5888 

295245 


397568  6 

1594323 


(4.17) 


where  we  have  set  q  =  4  and  used  (4.4)  to  simplify  the 
second  term. 

We  are  interested  in  finding  the  singularity  in 
the  Cy  curve  so  we  must  find  the  series  for  the  free 
energy  per  site,,  E,  and  the  specific  heat,  Cy.  Using  the 
Boltzmann  average  one  finds 


E  =  kt2  Ht  (ir ln  z)  ' 


and  of  course  we  know 


CTT  = 


dE 


V  dT 


(4.18) 


(4.19) 


It  is  desirable  to  express  these  quantities  as  derivatives 
with  respect  to  the  variable  x  and  w.  Using  the  definitions 
of  these  variables 


-J„/KT 


and  w 


1 

2 


J„/KT 
(e  H 


y 


X 


e 


(4.20) 


* 
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we  get 


d  _  JH  d 
dT  -  kt2  x  dx 

d  JH  /  ,  In  d 

dT  -  '  kt2  (w  +  2]  dw- 


(4.21) 

(4.22) 


Using  (4.15)  and  (4.21)  with  (4.18)  and  (4.19)  gives  the 
following  L.T.  series: 


E  =  2  JH  [4x9  +  6x8  +  21x{  -  12x8  +  54x9  +  150x10  -  242x1X 

+  940x12  +  2 86x13  +  ...],  (4.23) 

Cy  =  2K  (in  x)2  [l6x9  +  36x^  +  l47x ^  -  96x8  +  486x9 


+  1500x10  -  2662X11  +  11280x12  +  37l8x13  +  .  .  . ]  # 

(4.24) 

Equations  (4.17)*  (4.18),  (4.19)*  and  (4.22)  combine  to  give 
the  H.T.  series: 


F  _  o  t  rl  52^  200  2  241 6.  3  ,  13504  4  225280  5  ,  fll  p( 

E  -  2  JH  9  Hr*  729^  +  6561  +  5904 9W  531441  4  J 


Cy  =  4K  [ln(l+2w)  ]2[gy-  +  2l8'7W 


12  .  2 
656lw 


46120  3  ,  38528  4  n 

59049  +  531441  T  * ‘ *  * 

(4 . 2( 


We  note  that  due  to  the  nature  of  the  H.T.  variable  w  we  lose 
two  terms  in  going  from  |  In  Z  to  Cy  in  these  latter  series. 
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4 . 5  Estimation  of  the  Critical  Temperature  -  The  Pade  Approximant 

The  Pade  (47)  approximant  is  a  rational  approxima¬ 
tion  to  a  function  which  is  represented  by  a  power  series. 

To  a  given  series  f(z)  =  2  c  there  corresponds  a 

n  n 

unique  rational  fraction  F^  ^(z),  whose  numerator  is  a 
polynomial  in  Z  of  order  M*  and  whose  denominator  is  a 
polynomial  of  order  N.  When  expanded  in  powers  of  Z9 
F^  jyj(z)  agrees  term  by  term  with  the  series  f  (  Z)  for  more 
terms  than  any  other  rational  fraction  whose  numerator  and 
denominator  are  of  degrees  M  and  N  respectively.  The 
fraction  F^  ^(Z),  or  as  it  is  sometimes  written  [N*M]  *  is 
known  as  the  Pade  approximant  to  f (  Z) .  For  a  proof  that 
the  Pade  approximant  is  unique  the  reader  is  referred  to  a 
book  on  the  subject  of  continued  fractions  such  as  Wall  (44). 


The  approximant  is  of  the  form 


r  i  _  /„»  _  FyyWj  =  a°  +  aiz  +  a2 

[n,m]  fNjM(z)  bn^m(z) 


£  +  .  .  .  +  afIZM 


2  M 

bQ  +  b-^ Z  +  b^Z  +  .  .  .  +  a^Z  . 


(4.27) 


The  polynomials  AN  M(z)  and  BN  M(z)  are  determined  such  that 


f(z)  bNjM(z)  -  aNjM(z)  =  (zN  +  M  +  b  , 


(4.28) 


where  (ZN  +  M  +  1)  represents  a  polynomial  in  Z  whose  lowest 
term  is  ZN  +  M  +  1 .  one  can  sat;isfy  (4.28)  by  taking  the 


product 


57. 


cobo  +  ^coDl  +  clbo)  Z  + 

( c0b2  +  Clbl  +  c2b0 )  Z2  +  .  .  (4.29) 

M-f  TT> 

and  putting  the  coefficients  of  the  terms  in  Z  ,  where 

r  =  1,  2,  .  .  .  N,  equal  to  zero.  A  set  of  N  equations 

in  the  N  +  1  unknowns  bQ,  b^,  .  .  .  ,  b^  results  and  from 

this  we  can  find  b^,  b^,  .  .  .  b^  in  terms  of  bQ.  ¥e  then 

set  the  first  M  +  1  terms  of  the  R.H.S.  of  (4.29)  equal 

to  jvj(Z).  The  coefficients  aQ,  a^,  .  .  .  a^  are  all 

determined  in  terms  of  b  .  The  one  unknown,  b  ,  cancels 

o  o 

out  of  the  numerator  and  denominator,  and  we  have  thus 
determined  the  Pade  approximant  F^  ^(Z). 

If  the  function  f(z)  has  a  pole  at  some  point  Zq} 

then  better  and  better  approximations  to  Z^  will  be  obtained 

by  finding  the  zeros  of  ^(z)  for  larger  and  larger  values 

of  N.  Baker  (45)  has  used  this  method  to  estimate  the 

critical  point  from  the  series  expansions  of  the  derivative 

of  the  magnetic  susceptibility  for  the  Ising  model  on 

various  two  and  three-dimensional  lattices.  The  results 

obtained  by  examining  the  poles  of  the  sequence  of  [N,N] 

approximants  converge  rapidly.  For  instance,  for  the  simple 

-2J/KT 

quadratic  lattice  Baker  gets  x  =  e  G=  0. 41421058 

from  the  approximant  [7*7] ^  while  x  is  known  exactly  to 

0 

be  /2  -  1  =  0.41421356  ... 
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We  will  now  apply  the  Pade  approximants  to  the 
L.T.  and  H.T.  series  for  C^.  and  see  if  we  can  obtain  an 
estimate  of  the  location  of  the  critical  temperature.  This 
will  be  the  temperature  at  which  the  pole  in  the  specific 
heat  curve  occurs.  In  table  3  below  we  have  tabulated  the 
denominator  ^(x),  and  its  real  zero  which  corresponds 

to  the  pole  in  C  ,  for  the  [N,N]  approximant  to  the  L.T. 

Cy  series. 


[NjN] 

Denominator  B^  ^(x)  of  [N,N] 

X 

c 

[1,1] 

X 

1 — 1 

,-Q 

0 

[2,2] 

b  (25,065  x2  -  5,292  x  +  1296) 

No  real 
pole 

[3,3] 

b  (l2.220504x3  -  0.097543x2  -  O.629173X 

0  -  1.000000) 

0.47647 

[4,4] 

bQ (83 . 337721x^  +  12.080791x3  -  4.655109x^ 

-  7.660339x  -  1.000000) 

0.47968 

Table  3.  Results  of  the  PadlT  Method  on  the  L.T.  Series 


Since  we  do  not  have  enough  terms  in  the  series  to 
calculate  a  higher  order  approximant  we  must  base  our  estimate 
of  x  on  the  results  above.  After  examining  the  trends  in 


c 

convergence  in  Baker’ 
estimating  the  point 
three  figures,  to  be 


s  results,  we  feel  quite  safe  in 
of  critical  adsorption,  correct  to 


-Jh/KT 


C 


X  = 


ro 


e 


0.479. 


59. 


The  H.T.  series  for  Cy  has  been  calculated  up 
to  only  the  fourth  power  in  w.  Hence  the  highest  order 
approximant  which  can  be  found  is  [2,2].  Both  [1,1-}  and 
[2,2]  have  been  calculated,  but  unfortunately  the  pole  for 
[ 1 , 1 ]  is  at  a  large  negative  value  of  w  (which  has  no 
physical  significance)  and  the  poles  for  [2,2]  are  complex. 
To  estimate  the  critical  point  from  the  H.T.  series  by  the 
Pade  method  we  would  have  to  calculate  at  least  four  more 
terms  in  the  series.  The  graph  counting  involved  in  such 
a  calculation  would  be  very  laborious.  It  has  not  been  done 
in  this  work. 


4.6  The  Ratio  Method 


The  ratio  method  of  estimating  the  singular  point 

of  a  function  represented  by  a  power  series  has  been  used  by 

Wakefield  (43),  Ballentine  (8),  and  many  others  to  estimate 

the  critical  point  of  various  Ising  model  problems.  If  a 

series  2  A  converges,  we  know  from  the  ratio  test  that 
n  n  A 

the  limit  as  n  w  of  y  must  be  less  than  one.  If  the 

jH.  t 

n-1 

series  is  a  power  series  in  u  we  can  find  the  radius  of  con¬ 
vergence,  which  is  characterized  by  the  singular  point  closest 
to  the  origin,  by  taking  the  ratios  of  successive  terms  and 
extrapolating  to  find  the  limit  as  n  -  °o.  We  have 


n 


a 

i  •  i  n 
lim  I  — 

_>  oo  an-l 


1 

u 


(4.30) 


is 


where  is  the  coefficient  of  un  in  A  and  u 

n  n  c 

the  critical  point.  The  series  converges  for  |  u  |  <  u_  . 

If  we  plot  — H  vs.  p-  and  extend  the  graph  to  n  =  oo  9 

an-l 

we  can  estimate  u  ,  from  the  intercept  of  the 

c 

graph  using  (4.30).  In  this  method,  one  can  also  make  use 

of  the  fact  that  if  a  function  f(u)  has  a  pole  at  some 

critical  point  u  ,  then  successively  higher  powers  of 

c 

this  function  will  have  poles  at  this  same  point. 

The  L.T.  C  series  was  such  that  the  signs  of 
the  terms  were  mixed  positive  and  negative  but  were  not 
alternating.  Higher  powers  of  the  Cy  series  were  taken 
until  all  terms  became  positive.  However,  when  the  ratios 
were  plotted  on  the  graph  the  points  were  far  too  irregular 
to  make  any  extrapolation.  The  H.T.  Cy  series  had  so  few 
terms  that  only  three  ratios  could  be  found.  It  was  felt 
that  no  reliable  extrapolation  could  be  made  from  so  few 
points.  The  ratio  method  proved  to  be  unsatisfactory  for 
our  series  with  the  number  of  terms  presently  available. 


Chapter  5»  Discussion  and  Conclusion 

5. 1  Criticism  of  the  Various  Adsorption  Theories 

There  are  several  criticisms  that  can  be  made  of 
the  various  adsorption  theories  which  we  have  looked  at  in 
Chapters  II  to  IV.  Apart  from  the  other  simplifying  assump¬ 
tions  made,  the  monolayer  theories  may  be  criticized  simply 
on  the  grounds  that  there  are  so  few  systems  which  exhibit 
the  adsorption  of  only  one  layer.  In  the  Langmuir  theory 
all  lateral  interactions  are  neglected,  when  in  fact  we  know 
that  a  molecule  is  affected  by  the  presence  or  absence  of 
molecules  around  it.  The  treatments  of  Fowler,  Peierls,  and 
Ballentine  and  Betts  all  include  nearest  neighbour  lateral 
interactions  and,  although  it  can  be  shown  that  the  effect 
of  molecules  beyond  the  nearest  neighbours  is  small,  there 
is  bound  to  be  some  contribution  from  them.  The  first  two 
of  these  theories  use  approximations  (due  to  Bragg  and 
Williams,  and  Bethe  respectively)  and  are  known  to  be  in 
error  in  the  critical  region.  Because  the  treatment  due 
to  Ballentine  and  Betts  employs  only  exact  methods,  theirs 
is  certainly  the  best  localized  monolayer  theory  to  date. 
Some  indication  of  the  error  involved  in  using  the  approxi¬ 
mation  is  given  in  table  4,  which  compares  the  location  of 
the  critical  point  given  by  the  two  approximate  methods  with 


the  exact  location. 
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Method 

-J  /KT 

x  =  e 
c 

Bragg -Williams  Approx. 

0.3679 

Bethe  Approx. 

0.2501 

Exact 

0.1715 

Table  4.  Location  of  Critical  Point  as  given  by 
approximate  and  exact  methods  for  the 
Ising  Model. 

CernuschiTs  theory  for  two  layer  adsorption  ignores 
all  lateral  interactions ,  while  Dube’s  treats  them  in 
an  approximate  manner.  However,,  both  treatments  implicitly 
assume  that  there  is  an  interaction  between  the  gas  molecules 
in  different  layers,  since  the  condition  for  critical 
adsorption  is  that  the  second  layer  molecules  be  more 
tightly  bound  than  the  first  layer  ones.  A  strong  point 
of  the  Cernuschi  and  Dube  theory  is  that  the  sites  for 
molecules  adsorbed  into  the  second  layer  are  formed  by  four 
molecules  on  the  corners  of  a  unit  cell  (square)  in  the  first 
layer.  This  forms  a  close-packed  structure  which  is  the 
most  logical  way  for  adsorption  to  proceed. 

The  two  fundamental  objections  to  the  BET  theory  are 
firstly  that  lateral  interactions  are  neglected,  and  secondly 
that  the  molecules  are  adsorbed  one  on  top  of  another,  rather 
than  in  a  close-packed  structure.  Hill’s  treatment  was  an 
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attempt  to  remedy  the  first  point,  although  to  do  it  he,  too, 
used  the  Bragg -Williams  approximation.  BallentineTs  many 
layer  theory  is  good  in  that  its  treatment  of  lateral  inter¬ 
actions  is  mathematically  exact,  but  physically  it  only 
treats  nearest  neighbour  contributions.  However  this  theory 
is  based  on  a  rather  unrealistic  picture  of  adsorption  which 
permits  molecules  to  be  adsorbed  on  top  of  holes  in  lower 
layers . 

5. 2  Summary 

Several  theories  of  adsorption,  ranging  from  very 
elementary  treatments  of  monolayer  adsorption  to  more 
sophisticated  treatments  of  multilayer  adsorption,  have  been 
reviewed.  The  emphasis  has  been  on  localized  adsorption 
although  the  mobile  case  has  been  considered.  The  simplifi¬ 
cations  and  assumptions  made  in  each  theory  have  been  discussed 
and  the  weaknesses  of  each  have  been  criticized. 

A  model  for  localized  two  layer  adsorption,  in 
which  the  occupation  of  an  adsorption  site  is  considered  to 
be  one  of  three  possible  states  for  that  site,  has  been 
presented.  For  this  model  we  have  had  to  take  a  special 
case  to  facilitate  the  expansion  of  the  partition  function 
of  the  system  in  both  low  and  high  temperature  variables. 

The  effect  of  this  specialization  has  been  to  reduce  the 
number  of  free  interaction  parameters  from  three  to  one. 

Series  expansions  for  the  partition  function  have  been  found 


64. 


and  the  coefficients  have  been  calculated  using  the  graph 
counting  techniques  of  Domb  and  Sykes.  The  expansions  of  E, 
the  free  energy  per  site,  and  Cy,  the  configurational  specific 
heat  have  been  found  from  the  partition  function. 

Two  methods  have  been  used  in  an  attempt  to  estimate 
the  location  of  the  critical  point  of  the  two-layer  adsorbed 
system.  Only  one  method  was  successful.  The  specific  heat 
curve  has  a  pole  at  the  critical  point  which  will  show  up 
in  the  Pade  approximant  to  Cy.  The  critical  point  is  a 
singular  point  for  both  Cy  and  E,  so  these  two  series  will 
converge  for  |  x  |  <  x  .  In  theory  xc  could  be  estimated 
from  these  series  using  the  ratio  method.  In  fact  this  was 
not  the  case.  For  the  H.T.  series  too  few  terms  were  known, 
and  for  the  L.T.  series  the  ratios  were  far  too  irregular 
to  make  any  extrapolation. 

An  insufficient  number  of  terms  was  known  for  the 
H.T.  Cy  series  to  make  use  of  the  Pade  approximant.  However, 
for  the  L.T.  Cy  expansion  approximants  up  to  the  fourth  order 
could  be  calculated.  [3,3]  indicated  a  pole  at  x  =  0.4765, 
while  [4,4]  indicated  a  pole  at  x  =  0.4797*  On  the  basis 
of  these  results  we  feel  safe  in  estimating  x  only  to 
three  figures.  The  estimate  of  the  critical  point  of  the 
two-layer  adsorbed  system  which  we  give  is  xq  =  0.479  *** 
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To  compare  this  result  with  the  results  obtained 
by  Ballentine  for  the  application  of  the  Ising  model  to 
critical  adsorption,  we  must  note  that  the  variable  x  as 
used  by  Ballentine  is  the  square  root  of  the  x  used  in 
this  thesis.  Table  5  presents  the  locations  of  the  critical 
point  for  these  models. 


Model 

Location  of  critical  point  (x  ) 

Lx 

Ising  -  1  layer 

x  =  0.1716* 
c 

Ising  -  2  layers 

x  =  0.289 

0 

Ising  -  oo  layers 

X  =  0.4119** 
c 

Present  -  2  layers 

xc  =  0.479 

Table  5.  Location  of  Critical  Point  for  Ising  and 
Present  Models. 

Calculated  by  Kramers  and  Wannier  (28) 

■ft--#- 

Calculated  by  Wakefield  (39) 

For  a  constant  interaction  energy,  the  critical 
temperature  as  predicted  by  the  present  model  is  much  higher 
than  that  predicted  by  the  2  layer  Ising  model.  We  conclude 
that  critical  adsorption  will  take  place  in  a  system  of  gas 
molecules  which  have  energies  much  higher  than  the  cutoff  one 
would  expect  from  the  Ising  model. 
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5. 3  Suggestions  for  Further  Work 

The  best  estimate  of  the  critical  point  for  this 
model  is  based  on  the  fourth  order  Pade"  approximant  which 
uses  the  first  nine  terms  of  the  low  temperature  Cy  expan¬ 
sion.  If  one  could  extend  this  series  by  an  additional  eight 
terms ,  the  eighth  order  approximant  could  be  used  and  the 
value  of  xc  found  from  this  would  be  correct  to  about  six 
figures.  To  do  this  would  involve  counting  the  closed 
configurations  (both  Joined  and  separated)  of  up  to  and 
including  20  bonds.  Considerable  effort  would  be  required 
but  it  is  certainly  feasible. 

It  is  hoped  that  one  could  find  an  expansion  of 
the  partition  function  which  would  not  require  the  simplifying 
restriction  which  has  been  placed  on  the  interaction  parameters. 
If  this  were  possible  the  theory  would  be  a  completely  general 
one  with  three  parameters.  We  know  from  equation  (3.17) 
that  the  pressure  depends  on  the  chemical  potential  \i  .  The 
restrictions  given  by  equation  (4.4)  fix  the  value  of  \i  and 
hence  that  of  p.  Thus  to  calculate  theoretical  isotherms 
for  this  model,  it  is  necessary  that  we  make  the  generaliza¬ 
tion  and  remove  the  restrictions  (4.4)  However,  such  an 
expansion  will  be  very  complicated  because  terms  in  and 


are  involved. 
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The  two  layer  theory  was  originally  considered  with 
the  view  that  it  would  ultimately  be  extended  to  treat  three 
or  more  layers,  and  possibly  generalized  to  treat  n  layers. 

The  extension  to  more  layers  would  involve  choosing  an 
"occupation  variable”  with  four  or  more  possible  states.  This 
problem  is  suggestive  of,  but  not  the  same  as,  the  Ising 
model  for  particles  of  higher  spin  values.  In  this  respect 
it  is  interesting  to  note  that  Betts  (43)  has  recently  developed 
a  technique  for  treating  nearest  neighbour  lattice  problems 
where  each  element  has  four  possible  states.  His  techniques 
might  well  be  applied  to  the  three  layer  adsorption  problem. 

We  conclude  by  noting  that  the  model  which  has  been 
proposed  in  this  thesis  and  which  we  have  suggested  might  be 
extended  to  include  more  layers,  is  an  improvement  over  the 
Ising  model  for  several  layers  as  far  as  the  structure  of 
the  adsorbed  phase  is  concerned.  However,  it  may  still  be 
criticized  on  the  grounds  that  one  would  expect  the  adsorbed 
phase  to  have  a  close-packed  lattice  structure.  In  this 
respect  the  model  cannot  be  improved  because  of  its  inherent 
two-dimensional  nature. 


68. 


BIBLIOGRAPHY 


1.  de  Saussure,  T.,  Gilbert's  Ann.  der  Physik,  47,  113,  (l8l4). 

2.  Knudsen,  M. ,  Ann.  Phys .  Lpz . ,  34,,  593,  (1911 ). 

Ann.  Phys.  Lpz . ,  44},  64l,  (1915). 


3.  Ising,  E. ,  Z.  Physik,  31,  253,  (1925). 

4.  Fowler,  R.  H. ,  Proc.  Camb.  Phil.  Soc.,  32,  144,  (1936). 

see  also  Fowler  and  Guggenheim,  Statistical 
Thermodynamics.,  Camb.  Univ.  Press,  London, 
(1949). 

5.  Peierls,  R. ,  Proc.  Camb.  Phil.  Soc.,  32,  471,  (1936). 


6.  Bumble,  S.,  and  J.  M.  Honig,  J.  Chem.  Phys.,  33, 

424,  (I960). 

7.  Hijmans,  J.,  and  J.  de  Boer,  Physica,  21,  471,  485,  499, 

(1955). 

8.  Ballentine,  L.  E.,  Thesis,  Univ.  of  Alberta,  Edmonton,  (1962). 

9.  Ballentine,  L.  E.,  and  D.  D.  Betts,  Physica,  28,  1077,  (1962). 

10.  Jura,  G. ,  E.  H.  Loeser,  P.  R.  Basford,  and  W.  D.  Harkins, 

J.  Chem.  Phys.,  14,  117,  (1946). 


11.  Fisher,  B.  B. ,  and  W.  G.  McMillan,  J.  Chem.  Phys.,  28, 

549,  (1958). 

12.  Machin,  W.  D.,  and  S.  Ross,  Proc.  Roy.  Soc.,  A265, 

455,  (1961). 

13.  Singleton,  J.  H.  and  G.  D.  Halsey,  Jr.,  J.  Phys.  Chem., 

58,  ion,  (1954). 


14.  Polley,  M.  H.,  W.  D.  Schaeffer,  and  W.  R.  Smith, 

J.  Phys.  Chem.,  57,  469,  (1953). 

15.  Langmuir,  I.,  J.  Am.  Chem.  Soc.,  40,  1361,  (1918). 

16.  Fowler,  R.  H. ,  Proc.  Camb.  Phil.  Soc.,  31 ,  260,  (1935). 

17.  Williams,  A.  M. ,  Proc.  Roy.  Soc.  Edinburgh,  38,  23,  (1918); 

39,  48,  (1919).  Proc.  Roy.  Soc.,  A9 6, 

2H7,  298,  (1919). 


69. 


18.  Henry,  D.  C.,  Phil.  Mag.,  (6),  44,  689,  (1922). 

19.  Volmer,  M. ,  Z.  Phys.  Chem. ,  115,  253,  (1925). 

20.  Bragg,  W.  L. ,  and  E.  J.  Williams,  Proc.  Roy.  Soc., 

A145,  699,  (1934). 

21.  Cernuschi,  F. ,  Proc.  Camb.  Phil.  Soc.,  34,  392,  (1938). 


22.  Dube,  G.  P.,  Proc.  Camb.  Phil.  Soc.,  34,  587,  (1938). 

23.  Polanyi ,  M. ,  Verh.  Dtsch.  Phys.  Ges.,  16,  1012,  (1914); 

18,  55,  (1916). 

24.  Brunauer,  S.,  P.  H.  Emmett,  and  E.  Teller,  J.  Am.  Chem. 

Soc.,  60,  309,  (1938).  See  also 
Brunauer,  S.,  The  Adsorption  of  Gases 
and  Vapors,  Princeton  University  Press, 
Princeton,  (1943). 

25.  Halsey,  G.  D.  Jr.,  J.  Chem.  Phys.,  l6,  931,  (1948). 

26.  Hill,  T.  L.,  J.  Chem.  Phys.,  15,  767,  (1947). 


27.  Halsey,  G.  D.  Jr.,  and  W.  M.  Champion,  J.  Phys.  Chem., 

57,  646,  (1953). 

28.  Young,  D.  M. ,  and  A.  D.  Crowell,  Physical  Adsorption  of 

Gases,  Butterworths,  London,  (1962). 
Chapters  4,  5,  and  7- 

29.  Kramers,  H.  A.,  and  G.  H.  Wannier,  Phys.  Rev.,  60, 

252,  (1941). 

30.  Domb,  C.,  Advances  in  Physics,  9,  149,  (i960). 

31.  Newell,  G.  F. ,  and  E.  W.  Montroll,  Rev.  Mod.  Phys., 

25,  353,  (1953). 

32.  Hill,  T.  L. ,  Statistical  Mechanics,  McGraw-Hill, 

New  York,  (1956). 

33.  Fisher,  M.  E.,  J.  Math.  Phys.,  4_,  278,  (1963). 

34.  Onsager,  L. ,  Phys.  Rev.,  65,  117,  (1944). 

35.  Kaufman,  B. ,  Phys.  Rev.,  7 6,  1232,  (1949). 


70. 


3 6.  Kac,  M. ,  and  J.  C.  Ward,  Phys.  Rev.,  88,  1332,  (1952). 

37.  Sherman,  S.,  J.  Math.  Phys.,  1_,  202,  (i960). 

38.  Yang,  C.  N. ,  Phys.  Rev.,  8_5,  808,  (1952). 

39.  Van  der  Waerden,  B.  L. ,  Z.  Physlk,  118,  473j  (1941 ). 

40.  Wakefield,  A.  J. ,  Proc.  Camb.  Phil.  Soc.,  47,  419.,  (1956). 

41.  Oguchl,  T.,  J.  Phys.  Soc.  Japan,  6,  31*  (1951). 

42.  Domb,  C.,  Proc.  Roy.  Soc.,  A199,  199^  (1949) . 

43.  Wakefield,  A.  J.,  Proc.  Camb.  Phil.  Soc.,  47.,  799*  (1951). 

44.  Wall,  H.  S.,  Continued  Fractions,  van  Nostrand,  New  York, 

(1948).  Chapter  XX. 

45.  Baker,  G.  A.  Jr.,  Phys.  Rev.,  124,  768,  (1961). 

46.  Betts,  D.  D.,  Can.  J.  Phys.,  42  ,1564,  (1964). 

47.  Pade,  H.,  Thesis,  Ann.  de  l!Ec.  Nor.,  (3)*  9*  1*  (1892). 


APPENDIX 


71. 


CALCULATION  ON 


THE 


SEPARATED  LA 


TTICE 


CONSTANTS  OF  SIX  BONDS 


As  a  check  of  the  calculation  two  equations  are  given 


for  each  constant  but  the  last  one.  A  check  for  this  one  is 

provided  by  the  fact  that  the  total  number  of  configurations 

of  six  bonds  must  be  i  *7'  . 

\  o  / 


A. 


B. 


C. 


D. 


E. 


F. 


(  cf  ,/)  =  (O'  )(/>  -  5(  rf)  -  2(cf )  -  (CF)  -  2(dA 

-  2(l7  )  -  16N2  -  120N 

(□",/)  a  8(0  ,/)  -  (CO)  *  l6N2  -  120N 

( □  , A )  =  (O  )(A  )  -  8(0)  -  2 (CT)  -  (cC)  =  6M2  -  48n 

(□  ,A )  =  |  [6(0,/)  -  (CA )  j  =  6?J2  -  48m 

(Xs,/)  =  (X>)(/)  -  5()C  )  -  <X*)  -  2(0)  -  zafo 

-  2(-)-()  =  24k2  -  220N- 

(Xs,/)  -  1 2  ( X  ,  / )  -  (XV)  =  24M2  -  220N 
(X»A)  -  (X>(A)  -  6(X>  -  -  OO  -  (-H> 

-  ( u" )  -  6ll2  -  62N 

<X,A)  =  |[6(X,/>  -  (XV>]  =  6k2  -  62k 

(A,A)  =  |  (/AKA)  -  (A)  -  12( X )  -  2(><)  -  2(X\) 

-  2(C>)  -  2(y~( )  =  3n2  -  84n 

(A, A)  =  |  [  2(A,A)  -  2(>\)  -  (X)]  =  8s2  -  8ttN 
(AA/,/)  =  (AA/)(/)  -  5(/W)  -  2(AAA)  -  2( AAA  ) 

-  !/vV)  -  6(0>  *  2(C0)  -  (c£) 
c  2o4N2  -  2622N 

(AA/,/)  =|[6(AA,/)  -  2(OV)  -  2(AAA)] 

=  284N2  -  2622N 


72. 


c  •  ^  y~\  •>  / )  “ 

m* 

»•* 

( y~( ?/)  - 

H.  (V\Z_ ,/) 


(W.,/) 

i.  <\A/,/)  = 

(vA/s/)  = 

J.  (/vL  j  A  ) 


(AA  >A  ) 
K.  (AA,A) 


(AA,A) 
L.  (A/,  AO 


(AOA/> 


( )  (/ )  -  7(  )  - 


o  f  JV-/  > 


-  (  [A  )  -  (VuL  ) 


bi: 


O  *\r 

l2o:\ 


o 

">  /f  TvT  i- 
"j  Ui> 


■5  OAT 

J  aOx* 


%  j_2(AA ,/)  -  (  vvl)J 
«  (W. )(/)  -  ?(W_)  -  (AAL)  -  4(yX) 

-  (\Aa)  -  3  CO  )  -  2(c7)  -  4( 6/>  -  (aA /) 

=  l84n2  -  1708N 

=  3(AZ_,/)  «  2  (tf ,/)  -  (A/VL)  =  1S4N2  -  170SN 
(VV)(/)  -  5(vv)  -  (aA/)  -  2(\Az  )  -  2(00 ) 

-  3(vv )  -  2(0)  -  2<c0)  =  200N2  -  1820N 
-|[6(AZ-,/)  -  (/vO  -  2{lA,/)1=  200N2  -  1820N 
=  (AL)(A)  -  4(AZ_)  -  9 (Xs)  ~  2(vb)  -  8(>-<) 

-  l  n/\A  )  ”  \  /A'AZ—  i  *  VW  >  —  -  n  A  O' )  “  V  aA/ , 

-  (cO)  -  40  )  -  2 (AO  -  3(A)  -  4(eA > 

-  2(  O4 )  -  Adi-  )  »  2l6N2  -  2216N 

=  |[6(/sA,/}  -  (aAx)  -  2(00)  -  (A) 


-  (/W_)j  =  21610  -  2216N 
(AAUA)  -  3(AA)  -  2(A/V)  -  2CAAA) 

-  MW.)  -  2(vU  -  (aA)  -  3(A)  -  3(A) 

-  (AAA)  -  (00)  -  6(A)  -  2(A)  -  2(A) 

-  4(  A)  -  2(A)  -  6(0)  =  300N2  -  3128m 

|[6(/V\,/>  -  2(AAA)  -  (aA/)  -  3(vC0_ 

3COM2  -  3123n 

4  f  (A/)(A/)  -  (A/)  -  2CAA)  -  2(AA/) 

c  L 

-  2(AAA)  -  2(AA)  -  6(0)  -  4(00) 

-  10(  CO )  -  4(X)  -  2(l0)  -  4(lS)  -  4(DT) 

-  4(C0)  -  2(aA/)  -  12(D)  -  6(0)1 


✓  ™2 


162N 

ir 


--723N 


i[6(A/,A)  -  2(/W\)  -  (/W) 


=  162N2  -  172 


TT 


M.  (A y,/N)  “  (A/)(/X) 


N. 


0. 


D 

*.  * 


Q. 


*r> 

±  v  * 


«u 

6  ( Xs )  - 

<  W_  ) 

«cs» 

(□"X  - 

N» 

*•* 

m  • 

h 

2N2  -  772?? 

(AAA) 

1 

2 

1 - \ 

ON 

> 

-  (AAA) 

-  (Xv 

(□,/,/) 

«*o 

U9* 

1 

•<<» 

A 

[<□,/)(/ 

)  -  ?(□ 

,/) 

-  (c/ 

(□"X]  = 

2?J3  -  31 

N2 

4-  126N 

<D,/,/) 

- 

(/,/)(□)  - 

2(D)  - 

(C 

X)  -  ( 

4(D,/)  - 

in' ,/) 

•»« 

2N3  -  3 

(AA,/, 

/) 

M 

A  i  (  /v\  j 

/)(/)  - 

?(AA,/) 

Y3. 

(AA)  -  (W-)  -  (AAA)  -  (X"7) 

_ Z.  .  ,  r-.A 


2( 


(M,/,/) 


(AA,/,/)  » 


-  2(W_,/)  -  (\A/,/)  "  2(M,A)  -  (aA/) 

-  2(AAA)  -  3(v0)  -  2(D/s/)j 
ICON3  -  19XSN~  +  96721? 

v;  i  6 ( A/ s/s/)  “  Sv  □  ,/,/)  -  2 { AA/ s / )  J 

100N3  -  1918N2  4-  967  2N 

|  [(AA ,/)(/)  -  5(/\A,/)  -  AAA.,/)  -  (W.) 

-  WX./)  -  3<X'  »/)  -  2(vA/s/)  -  2( AA , A ) 

-  2('XX)  -  (aA/)  -  (AAA)  -  2(CrV)l 


72K-5  -  1372N' 


6856N 


(AA,/,/) 
(>(X,/)  = 

(>{,/,/)  = 
(A/, A,/) 


=  i[4(/V,/,/)  -  2(vV,/)’  -  72N3  -  1372N2.>  68?6n 

|[(X,/)(/>  -  5(  X  >/)  -  ()<",/)  -  2(>(,A) 


(W)l  -  2K3  -  39N2  -6-  198k 


fr  [(A,/,/)  -  <)\X/)j  *  2N3  -  39lX  f  198N 
=  (A/,A)(/)  -  5(AAA)  -  2(AA,A)  -  (aA/) 

-  2(AA,A)  -  4(/V,/V)  -  3(A/,/f0  -  (\Az.) 

-  >r(a,A)  -  2(/W\)  -  (AAA) 

=  2l6N3  -  4324-K2  <■  22996N 

=  i  fl2(A/,/,/)  -  ^(A/ s A/ )  -  2(A/V,/) 

2  L 

-  2(v\/,/)  |  =  2l6w3  -  I+324U2  +  22996N 


(/V,A  ,/) 


s.  (A,  A,/) 


v . 


T/F 


(A, A  ,/) 


Y4 . 

(AN,  A  -fv/';  -  5(A\,A  )  -  4(/\,A  .  -  ( AZ_  j A  ) 

-  2(/,\,A/)  -  6(/f\,/N)  -  (>Cy)  -  (/W_) 

-  (>-<)  -  2()/X)  =  48?;3  -  968s2  +  5168N 

i  [i2(A,/,/)  -  2(A,A/)  -  (W.,/) 

-  ()<y  >/)  j  =  48n3  -  968n2  +  5168N 
<A,A,A)  =4  r(A,A)(A)  -  2(A,A>  -  2(A/,A)  -  (X) 

-  2(A/\,A)  -  6(A) A )  -  6(X)A>  -  2(AL,4) 

-  4(0, A)  -  2(AA/)  -  (AAA)  -  3(W_) 

-  (AAA)  -  (aA/)  -  3(X^)  "  3<X>  '  (C^) 

-  4(>-()  -  (o4)  -  (XX)] 

—  O  <7',  )  ,  >r ^  ) ,  -|  } ,  ).  M 

—  J)0-j  —  /,r;J  *J>  irJ.tTxriv 

(A,  A,  A)  -K6<A,A,/)  -  2(AA,A)  -  ( AZ_  ,  A  >  j 


o  u 


36s3  -  744nc  -v  4l44N 


2 


U.  (A/,/,/,/) 


1  r 

“zr  * 

J  L 


(A/, /,/)(/)  -  5( A/,/,/)  -  2(M,/,/) 


-  2(AA/ ,/)  -  2(A/AO  -  2(A/SA,/) 

-  2(AA>/>/)  ~  2(vW>  “  4(  □,/,/)] 

24k4  -  71 2N3  +  734 2N2  -  26604k 

1  r 


(A/,/,/,/)  =  3  l6(A,/’/’/)  "  2(A/\,/,/)j 

=  24rr  -  712A  *  7342N2  -  26604N 

(/..,/./,/)  =  ••  [(A, /,/)(/)  -  5(A,/,/)  -  4(X,/,/) 


2(A; A  ,/)  -  (A, A/)  -  (W.  ,/)  -  C>(\/) 


,  ,  ,  ,0  _  1'  4  w„.3  42.88„2 

-(AL)/)/)  j  “  -  loOrX  v  — a — N  -  6 

•1 


6o48n 


•?  r2(A, /,/,/)  -  (AA  ,/,/)] 
loON3  +  — A 

*  A 


(A, /,/,/) 

1A4 

c,  J.  A /  «T  I 

J 

(A  ,A  ,/,/)  =  c  [i8(/\, /,/,/)  -  2(M,/,/)  -  (AA,/,/) 


A 

o 


04- Bit 


2(/V,A,/)j-  36 


IT'-  lO^N3*  11600N2-  43446N 


X. 


Y 


(A  ,A  ,/,/)  -  :2  ;_(A, A ,/)</)  -  ?(A,A,/)  -  2</V,A ,/) 

-  3(/1\,A  ,/)  -  6(A,A,A)  -  (A/V,/) 

-  2(/V\,A)  -  (V\X_,/)  -  (X>/)  -(AA,A)] 

=  3oNv  -  I094N3  +  lloOON2  -  43446M 

(A,/,/,/,/)  r[(  A, /,/,/)</)  -  5(A, /,/,/)  -  3(/K»/,/,/) 

-  2(/V, /,/,/)  -  4(A,A,/,/)  -  2(/v\,/,/) 

-  (A JL,/,n  -  (A/,A,/)J 

=  4’!'’ 

(A,/,/,/,/)  =  |  [30(/, /,/,/,/)  -  2(A/,/,/,/)] 

=  4N''  -  164N  ’  2607n3  -  X9217N2  --  56040tf 

(/,/»/»/,/,/)  *  |  [(/,/,/,/,/)(/)  -  5(/, /,/,/,/) 

-  2(A ,/,/,/,/)  -  (/V, /,/,/)] 


l64li  ’  *  25 07 N 3  -  19217N2  *  56040N 


_  4,, 6  14..-5.  905-4  2239-3 .  584627-2  47644- 

-  irj’A  -  -j-.  *  ~  — 2“ ^ — a - t  ^ 


Pinal  Chock: 

Total  connected  constants  of  six  bonds  -  1628N 


Total  separated  constants  of  six  bonds 
b„6  2„5  ,  17,A  5.T3  .  137^2 

XT’. .  -  +  -Q-; J  ~  r:i  -qq-iM 


5! 


v 


^-885. 


T\T 

■N 


Total  number  of  configurations  of  six  bonds 


